Draft version February 4, 2008 

Preprint typeset using I4TgX style emulateapj v. 6/22/04 



ROCKING THE LIGHTHOUSE: CIRCUMPULSAR ASTEROIDS AND RADIO INTERMITTENCY 

J. M. Cordes & R. M. Shannon 

Astronomy Department, Cornell University, Ithaca, NY 14853 
Draft version February 4, 2008 

ABSTRACT 

We propose that neutral, circumpulsar debris entering the light cylinder can account for many time-dependent 
pulsar phenomena that are otherwise difficult to explain. Neutral material avoids propeller ejection and injects 
sufficient charges — after heating, evaporation, and ionization — to alter current flows and pair-production 
and thus trigger, detune, or extinguish coherent emission. Relevant phenomena, with time scales from seconds 
to months, include nulls, rotating radio transients (RRATs), rapid changes in pulse profile ("mode changes"), 
variable subpulse drift rates, quasi-periodic bursts from B1931+24, and torque variations. Over the ~ 10 7 yr 
lifetime of a canonical pulsar with 10 12 G surface magnetic field, less than 10 22 Bi2 g (< 1O~ 6 M0) is needed to 
modulate the Goldreich- Julian current by 100%. Circumpulsar material originates from metal-rich, supernova 
fallback gas that aggregates into asteroids. Debris disks can inject sufficient material on time scales of interest, 
yet be too tenuous to form large planets detectable in pulse timing data. Asteroid migration results from 
collisions and the radiation-driven Yarkovsky and Poynting-Robertson effects. For B1931+24, an asteroid 
in a ~ 40 day elliptical orbit pollutes the magnetosphere stochastically through collisions with other debris. 
Injection is less likely for hot, young and highly magnetized pulsars or millisecond pulsars that pre-ionize any 
debris material well outside their small magnetospheres. Injection effects will therefore be most prominent 
in long-period, cooler pulsars, consistent with the distribution of relevant objects in P and P. A pulsar's spin 
history and its radiation-beam orientation may influence whether it displays nulling, RRATs and other effects. 
Subject headings: stars: neutron — pulsars: general — pulsars:individual (B 193 1+24, B0656+14) — accretion 
— acceleration of particles 



1. INTRODUCTION 

Rotation-driven pulsars display myriad phenomena in their 
radio emission on time scales that range from nanoseconds 
to years. Here we consider a unifying model for phenom- 
ena that occur on time scales from pulse-to-pulse to month- 
like or longer, including: (1) null pulses, where the emis- 
sion evidently vanishes but returns in pulses that maintain 
synchronism with the spin; (2) drifting subpulses, whose 
systematic variations in pulse phase appear to be related to 
circulation of emission regions around the magnetic pole; 
(3) mode changes, where the pulse shape changes discon- 
tinuously between two or more shapes, often accompanied 
by dramatic changes in flux density; and (4) rotating ra- 
dio transients (RRATs) that correspond to bursts of single 
or a fe w pulses spaced by very long intervals (minutes to 
hours) ( McLaughlin et al. 2006). A few objects show giant 
pulses that can b e > 10 3 larger in amplitude than the av- 
erage pulse (e.g. ICognard et alj H9961 iHankins et al ] 120031 
Johns ton & Romanil2003UCordes et al.l2004HCairnsl 2004 ) . 

Coherent emission wanes on time scales ~ 10 7 yr for 
canonical pulsars — those with surface magnetic fields ~ 
10 12 G — and on much longer times for millisecond pulsars 
(MSPs; 10 s — 10 9 G). We assume, of course, that the spin- 
down time scale, r, = P/2P, is a reasonable, but by no means 
perfect, proxy for chronological age. The distribution of pul- 
sars in the period-period derivative (P-P) diagram shown in 
Figure [2 displays a dearth of pulsars with P < 10~ 16,7 ss~ l P 3 
in the bottom right corner that reflects these time scales. 
However, it is not clear whether pulsars shut off their ra- 
dio emission according to a threshold effect, such as the in- 
ability to sustain electron-positron cascades whose flows in- 
duce coherent emission, or that the rotational energy loss rate, 
E = /OO, simply cannot provide sufficient power to drive 



the ra dio beam luminosity jArzoumanian. Chernoff & Cordesl 
120021) . Empirically, the latter seems favored because the den- 
sity of pulsars in the P-P diagram rolls off smoothly rather 
than sharply at a "cliff". However, pair cascades may ter- 
minate differently in pulsars with the same P and P if sur- 
face magnetic fields show object-dependent complexity, thus 
smearing out any cliff effect. 

In this paper we develop a mechanism that pertains to the 
four phenomena described above but is particularly focused 
on nulling pulsars and RRATs. Both show episodes of ra- 
dio quiesence that appear to be totally devoid of any radio 
emission to levels < 10~ 2 of their 'on' states. Off states can 
last for long time periods (from minutes to days), although 
some nulling pulsars are off only for a few pulse periods. Both 
classes of objects are largely found on the right-hand periph- 
ery of the population of canonical pulsars in t he P-P dia - 
gram (Figure^. As pointed out by many (e.g. lBiggslll992l) . 
nulling may signify a faltering emission process that signi- 
fies incipient complete shut off of radio emission. Instead, 
we suggest that nulling and RRATs are both caused by pol- 
lution of NS magnetospheres by low-level accretion from a 
circumpulsar asteroid disk. This idea is motivated by the fact 
that one nulling pulsar, B 193 1+24, has state durati ons of days 
that a re cyclical with a ~ 40-day q uasi-p eriod ( Kram eret alJ 
2006). In addition, early work by Cheng ( 1985) explored the 
role played by interstellar grains in the current flows of active 
pulsars and their manifestation on the spindown rates of pul- 
sars. Our model differs by considering material that originates 
in a reservoir of orbiting debris, most likely supernova fall- 
back material that has formed a low-mass asteroid belt. Or- 
bits around neutron stars provide a plethora of time scales that 
can be imprinted on electrodynami c processes by material in- 
jected into the magnetosphere. Michel & Dessler ( 1981) con- 
sidered plasma disks around pulsars to address issues of elec- 



2 



Cordes & Shannon 



trodynamics in pulsar magnetospheres. Here too, our model 
is different in considering a largely neutral disk comprising 
macroscopic objects. 
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FIG. 1. — The period-period-derivative diagram for pulsars. Rotating 
radio transients (RRATs), pulsars that show giant pulses, and nulling pulsars 
that display different nulling fractions are distinctly labelled (see legend). 
Points labelled "UL" indicate pulsars for which upper limits have been placed 
on the nulling fraction. Solid lines of constant magnetic field and lines of 
constant characteristic spindown age, P/2P are shown. The most leftward 
dashed line denotes constant magnetic field strength at the light cylinder (with 
a field strength equal to that of the Crab pulsar) that seems to be a trademark 
of objects showing giant pulses (e.g. Cognard et al. 1996). The inner-gap 
"death line" is drawn as the rightward dashed line at an arbitrary location 
that is consistent with the diminishment of the radio pulsar densi ty in the 
diagra m. The middle dashed line is an outer-gap death line from Ta kata et alj 
1 2006). The pulsar B0656+14 is circled because it is an object with high 
surface temperature that shows low-level nulling. 



Cheng's (1985) grain inflow model was developed to ex- 
plain an apparent correlation between PP and space veloc- 
ity that was identified in the 1980s. While that correlation is 
now thought to be insignificant for canonical pulsars 1 , extrin- 
sic charges from inflowing grains also solved the issue of how 
a NS avoids being charged up by self-generated current flows 
and it also obviated the need for self-sustained pair produc- 
tion. However, the grain model appears to work only for cool 
(~ 10 s K), long-period pulsars because grains are evaporated 
by thermal radiation at distances of 10 9 - 10 10 cm, compara- 
ble to the light-cylinder radii of long-period pulsars, but much 
larger than those of short-period and, presumbably, younger 
and hotter pulsars. For our purposes, the longer-period and 
(presumably) cooler NS are just the sample we wish to ad- 
dress. In many respects, Cheng's analysis also applies to 
grains delivered from asteroids that are evaporated and ion- 
ized inside the magnetosphere. The differences are that (a) 

1 Note that space velocities are much less for pulsars that have under- 
gone accretion-induced spinup and consequent diminishment of the dipole 
field component. However, there does not appear to be a correlation between 
Z?dipole oc PP and space velocity among pulsars that have not accreted much 
material. 



orbital delivery is expected to be episodic, as opposed to the 
steady flow of interstellar grains; (b) asteroids and rocks will 
penetrate further into the magnetosphere than the ^0.1 /im 
grains considered by Cheng; and (c) circumpulsar disks prob- 
ably vary widely from pulsar to pulsar, providing a natural 
explanation for why pulsars with very similar spin parameters 
(P, P) are dramatically different with respect to the occurrence 
of nulls, etc. 

In the next sections, we describe relevant phenomena and 
key elements of the model for producing nulls, RRATs and 
other bursts, including (a) formation of the circumpulsar disk; 

(b) delivery of asteroids from the disk to the magnetosphere; 

(c) tidal disruption, evaporation, and ionization; and (d) ef- 
fects on gap accelerators and radio emission. We also discuss 
implications for spin variations of NS from the torque fluctu- 
ations that are expected and we suggest possible tests for the 
model. 

2. SALIENT PROPERTIES OF NULLS AND RRATS 

The most striking feature of both nulls and RRATs is 
the rapid transition between states, either on to off or vice 
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thus too fast to be caused by propagation effects in the 
interstellar medium (ISM). The nulling phenomenon includes 
pulsars that are mostly on and occasionally show nulls with 
durations of a few spin periods and others where the pulsar 
is mostly in the off state for periods of up to days. In a 
few well studied objects, a slow roll-of f of the pulsed flu x 
precedes a null. In one case (B 1944+1 7: iDeich et aT1(ll986ly ) 
the slow decay is over a few pulse peri ods and in another, 
a few tens of periods (J1752+2359; Lewandows ki et all 
(2004)). Nulling objects occur throughout the grouping of 
canonical pulsars in the P—P diagram. However, those with 
null fractions > 20% are strikingly more common near the 
right-hand side of the grouping. Notable outliers exist, such 
as B0656+14, as labelled in F igure ^ Data used to form 
Figure 1 are given in Table lAll and are from multiple sources 
in the literature and also from the ATNF/Jodrell pulsar catalog 
Qhttp : / / www .atnf.csiro.au/ research/pulsar /psr cat 7\ 

Some objects with nearly the same period and period 
derivative null while others do not. This suggests that P and P 
are not the sole determinants of a pulsar's instantaneous radio 
luminosity. The surface magnetic field strength is commonly 
estimated to be B n = (PP-i 5 ) l/2 (with P = 10" 15 P_i 5 ) and 
numerous work has associated either the pseudo-luminosity 
SD 2 , where S is the flux density averaged over period and 
D is the distance in kpc, or the true b eam luminosity (e.g. 
lArzoumanian. Chernoff & Cordesl 120021) with a power-law 
model P x P y , with exponents x,y determined from some fitting 
process. The presence and absence of nulling objects with 
the same spin parameters suggests that other variables are in- 
volved. Intrinsic models might include instabilities in the cur- 
rent flow th at are regulated by the thermal time con stant of the 
NS crust ( Cheng & Rudermanl ll980t lJonesl[l982l) or incom- 
plete filling of the open field line region with emission r egions 
associated with drifting subpulses (e.g. Deshnand e & Rankinl 
1999). 

Figure |2] shows the nulling fraction plotted against the 
angle between the spin axis and the magnetic axis a for 
the known nulling pulsars. Spin-magnetic axis angles 
were compiled from the analyses of Rankin (1990 , 1993) , 
lEverett & Weisberd J200l . and iMitra & Rankinl J2002I) . 
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With the exception of B 2315+ 21, which has a nulling frac- 
tion of 3.0 ± 0.5%(Biggs 1992) , all the nulling pulsars have 
a spin axis-magnetic axis offset of at most 60°. 
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FIG. 2. — Pulsar nulling fraction as a function of inclination angle a be- 
tween the magnetic and spin axes. The model presented in this paper suggests 
that orthogonally rotating (a ss 90° ) would be less likely to null, as the beam 
would severely disrupt asteroid transport into the magnetosphere, which is 
consistent with these observations. 



RRATs are perhaps an extreme form of nulling. Pulse event 
rat es range from 0.3 to 20 h r" 1 in the eleven objects reported 
bvlMcLaughlin et alJ J2006I) . A twelfth object dCordes et alJ 
2006) has a much higher rate (~ 200 hr" 1 ) was also discovered 
as an aperiodic source but later found to be periodic. 

RRAT pulses typically occur singly, in contrast to nulling 
pulsars in their on states. Though RRATs are biased to long 
periods, as with nulls there is a short-period example. The 
RRAT Galactic population is probably very large in spite of 
the small, known sample because their low duty cycles im- 
ply a small probability of detection and because the Parkes 
multibeam survey covered only a narrow range of Galactic 
latitude. Non-pulsed X- ray emission has been detected from 
the RRAT J1819- 1458 (iRevnolds et all2006l) . which has the 
largest apparent magnetic field (5 x 10 13 G), and has the low- 
est spindown age (117 kyr) of the 1 1 RRATs. 

3. BASIC PICTURE 

The essential feature of our model is a disk of asteroids 
that is sustained in largely neutral form for tens of millions 
of years from which asteroids or their fragments are injected 
into the magnetosphere of a NS. To account for observed nulls 
and RRATs, injection must be episodic but, over long time pe- 
riods, statistically ho mogeneous. Thu s, our model differs sub- 
stantially from that of lChengl i ll 9851) . who considered steady 
injection of grains from the general ISM, requiring their sur- 
vival through the bow shocks that we know are associated 
with essentially all pulsars, owing to their large space veloc- 
ities, Vp Sr ~ 10 2 - 10 3 km s" 1 . Moreover, his model implies 



that pulsar radio emission should always be steady (or al- 
ways quenched) rather than intermittent. Cheng's model may 
be relevant for older, cooler neutron stars that have no de- 
bris disks around th em. Our model also differs substantially 
from that of Colgate & Petschek ( 1981), who considered in- 
terstellar objects that would reach the NS surface in producing 
gamma-ray bursts from Galactic NS. 

A supernova fallback disk (see below) is metal rich and can 
be sufficiently low mass (1O~ 6 M ) or compact (inside the tidal 
disruption radius) to prevent formation of massive planets that 
would be detectable individually in pulsar timing measure- 
ments. However, rocks can form that will resist tidal shred- 
ding until they are well inside the magnetospheres of long- 
period pulsars and will stay neutral until they are melted and 
ionized by radiation from the NS surface or magnetosphere. 
Delivery of rocks to the magnetosphere can occur from colli- 
sions and orbital perturbations by protoplanets if the disk ex- 
tends beyond the tidal disruption radius. Alternatively, rock 
migration from a compact disk will occur from the radiation 
driven Yarkovsky effect. 

The situation we identify is essentially the complement of 
that needed to account for planet formation around ordinary 
stars. Planet formation requires that planetesimals exist out- 
side the gravitational tidal radius r tg ~ 10 11 cm (see below) 
so that runaway accretion followed by oligarchic growth and 
mergers of oligarc hs can form planetary core masses o f Earth 
mass or larger (e.g. Rafikov 2003; Goldreich et al. 2004j). The 
planet building process takes only ~ 10 5 yr in the inner solar 
syste m but cleanup o f debris can take several hundred million 
years (Goldrei ch et all2004l) . For our case, the disk need not 
be stable; in fact it can disperse on a time scale ~ 10- 100 
Myr, the apparent radio-emitting lifetime of a canonical pul- 
sar. 

Once a disk becomes neutral, asteroid migration will follow 
from many of the processes that have been identified in the 
solar system. These include formation and growth of plan- 
etesimals through collisions, destructive collisions, tidal cir- 
cularization and spin alignment, shepherding by larger bod- 
ies, injection into resonances, heating from photon and parti- 
cle fluxes and subsequent anisotropic re-radiation, and heat- 
ing from induction and tidal flexure. Once injected into the 
magnetosphere of a pulsar, asteroids will evaporate and, when 
ionized, provide charges that are subject to strong electromag- 
netic forces. Such forces are most prominent in "gaps" where 
there is a potential drop along magnetic field lines that will ac- 
celerate charges inward and outward. Even though the struc- 
ture of such gaps has been a subject of analysis for more than 
three decades, a definitive picture has yet to emerge. How- 
ever, the prominence of pulsars at both radio frequencies and 
in high-energy bands attests to the fact that such gaps exist. 
In Figure|3]we show a schematic magnetosphere for an object 
with a spin-magnetic-moment angle a = 20°. The gaps shown 
are schematic representations that are consistent with those in 
the literature but may underestimate their true extent if gaps 
extend all the way to the inner gap. It is plausible for the ge- 
ometry shown that asteroidal material that penetrates the LC 
will be subject to gap potentials. 

Injected charges are accelerated in gap electric fields and 
can activate a quiescent region, such as an outer-gap, as a 
source of gamma-rays that pair produce, or they can temporar- 
ily shut off the pair production of an active gap by reducing 
the accelerating field. A pair avalanche from an outer-gap (as 
opposed to a sustained cascade) can infl uence o t her r egions 
in the magnetosphere, as suggested by Wright (2003)and 
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IZhang. Gil & Dvksl £2006) for self-driven processes within 
the magnetosphere rather than by extrinsic causes, as we ar- 
gue here. 

In order for this picture to be viable, it must satisfy three 
constraints: (1) the disk must be sufficiently low-mass that it 
not produce planets that would be detectable in pulse timing 
measurements or by other means; (2) the rate of mass inflow 
from the disk into the LC must be large enough to influence 
magnetospheric current flows on time scales of interest; and 
(3) the disk must survive long enough (e.g. > 10 Myr) to in- 
fluence the electrodynamics of old pulsars. Next we set limits 
on the properties of the disk with these constraints in mind. 

3.1. Observational Constraints on Disk Masses 

Infrared observations with S pitzer and previous instru- 
ments (e.g. Brvden et. al. 2006, and references therein) have 
failed to detect disk emission around the MSP B 1257+12, 
which harbors at least three planets llKonacki & Wolszczanl 
120031) . Nonetheless, the observations still allow a debris disk 
around the MSP consisting of (e.g.) 100-km asteroids with 
as much as 10~ 2 Mm total mass. Related the oretical work 
jMiller & Hamilton l200lllBrvden et. all2003) imply that as- 
teroids larger than ~ 1 km can exist for longer than 1 Gyr 
a round the MSP. 

Phillips (1993) placed limits on asteroidal masses around 
the Vela pulsar from an upper limit on the unpulsed flux S u 
that would be produced by reflections of the pulsed radio flux 
S p from an asteroid belt. His limit for a power-law size dis- 
tribution of slope -3.5 extending from R„2 =100 km down to 
R a [ equal to the 20 cm radio wavelength is 
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For compact disks with tau ~ 10~ 2 , this result becomes mod- 
estly constraining as M a < 10~ 2,3 M® but does not limit the 
mass in rocks that are smaller than the 20 cm wavelength. 

Pulsar timing measurements are sensitive to orbital 
perturbations from objects of about the Moon's mass 
( Konacki & Wolszczan 2003). The rms sinusoidal timing 
residual caused by a single object of mass m a in a circular 
orbit with period P or b is 
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Limits on sin gle planet masses de pend on data set lengths but 
are typically (Thorse ttet all 19931) ~ lM e forP orb ~ 1 yr. An 
ensemble of N a asteroids in a belt with total mass M a with rms 
semi-major axis a rms yields an rms timing residual 
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where cz rms is in AU and £ w 2 is the ratio of rms to mean 
asteroidal mass. 2 

Timing residuals for long-period pulsars are typically in the 
range of 0. 1 to 1 ms (rms) and show a variety of characteristic 
time scales as well as displaying non-stationary behavior (i.e. 
rms residual that depends on data span length). The limit on 
the total mass in an asteriod belt is therefore 



M a /M m < (0.1-1) X 
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2 For a pulsar timing data set of length T, only orbital periods < 2izT 
should be included in the calculation of <5z ims because a low-order polynomial 
that is fitted and removed filters out the contributions of longer-period objects. 



Thus disks of 3> IM® aggregate mass in N„ ^> 1 aster- 
oids can easily reside around many pulsars, particularly if 
firms "C 1 AU, as we consider here. In addition, such disks 
can contribute significantly to the timing noise while not 
yielding detection of any individual asteroid. To detect at 
the £ x a level the sinusoidal variation in a Fourier analy- 
sis of «</ = n^jlO 3 data points with white noise measure- 
ment errors of cr W) i ms, an asteroid would have to satisfy 

flATjsin i(m a /M®) > 0.04£a w ^n d -{ . In the case where timing 
residuals are dominated by asteroidal noise from N a asteroids, 
to be detectable an asteroid must have a mass that is a multiple 
of the mean asteroidal mass, (m a )M a /N a , of 



J ( m a) £ ^C(a r ms /flmax) V N a / n d > 



(5) 



which is large for a large number of asteroids and reasonable 
values of other parameters (e.g. £ = 2, a lms /a mRX = l/2> and 
nd i=s 10 3 ). We conclude that canonical pulsars can harbor 
asteroid disks with a sizable fraction of an Earth mass in the 
disk and show timing noise consistent with that measured. As 
shown below, very low mass asteroid disks are sufficient to 
disrupt the magnetospheres of these same pulsars. 

3.2. Magnetospheric Currents and Disk Lifetime 

To see if a disk can provide mass at interesting rates after 
10 Myr, we consider the current and implied mass rate needed 
to alter the electric field in a magnetosphere. In order for the 
acceleration region to exist, it is necessary for the magneto- 
sphere to have low density. 

The fiducial charge densit y in such objects is of order 
the Goldreich- Julian density (Goldreich & Julian 1969j), p e = 
— f2 -B/lirc, corresponding to a number density 
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Here Z is the charge per particle, ip is the local angle between 
the spin axis and dipolar field, and P* is the NS radius. 

The current from the magnetic polar cap (defined as the 
footprint of dipolar field lines that do not close within the LC) 
with area Ape is (in charges s" 1 ) 



Moj = c«gjApc : 
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For fi p proton masses per charge, the corresponding mass rate 
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only a fraction 10~ 116 of the Eddington mass accretion 
rate for a NS. Note that the entire magnetosphere has a 
minimum mass M mag ps 47rP sur f ace P 3 ^ / ,m ; ,ln(rLc/P*)/ecP ps 
UfgBnUpP- 1 . 

For a pulsar born with a 10 ms spin period that spins down 
to 10 s in 10 7 yr, the total mass that would have to be accreted 
to account for the GJ current is, if we assume constant mag- 
netic field, 
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which is equivalent to the mass in a ~ 100 km sized asteroid. 
dMichel & DessleJT9 8 1 ) made a similar point in their discus- 
sion of much more massive plasma disks around pulsars. A 
pulsar cannot accrete during approximately the first 10 5 yr of 
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FIG. 3. — The magnetosphere of a pulsar showing the locations of outer- 
gap acceleration regions where E\\ The asteroid belt is shown (shaded) 
in the z = plane. The other shaded regions show outer gaps whose extents 
along and transverse to the last closed field line within the light cylin der (LC) 
are arbitrary but consistent with discussion in the literature (e.g. Ta kata et all 
2006). These gaps may extend all the way from the LC (denoted as vertical 
dashed lines near the horizontal boundaries of the diagram) to the NS surface, 
and the inner gap (too small to be shown). The spin axis is the arrowed 
(vertical) dashed line and the magnetic axis is shown by a arrowed (slanted) 
solid line. Length scales are in units of the LC radius. 



its life because tidal shredding and thermal radiation will ion- 
ize material before it can enter the magnetosphere. However, 
a pulsar spends most of its lifetime at long periods, so our 
estimate is relevant. 

Consider a disk that contains a total mass M a in N a asteroids 
within minimum and maximum radii, r\ , T2 ^S> r\, and thick- 
ness 2H = e/,r2. Inward migration implies a lifetime (for fixed 
N a ) 

M a N a 

Td= 7T = lT- (10) 

M a N a 

The flux of asteroids through a cylinder centered on the NS 
with fiducial radius r and surface area A c = lire^r is N a = 
n a A c \r\ 2N s .r\r\r^ 1 , implying a disk lifetime that must be 
larger than the radio lifetime T psr . This yields the constraint 

^<r£ (ID 
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In order to provide at least a fraction e G j of the Goldreich- 
Julian mass rate, we require M a > e G jM G j, where M a = N a m a 
and m a is a typical asteroid mass, or 



2r\r\ e G jM G j 
ft M„ 



(12) 



Evaluating the two constraints assuming, while taking r psr : 
100 Myr), we find 
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where we have assumed uniform mass density p in an aster- 
oid. It is obvious that a disk can easily satisfy these constraints 
and not be detectable in pulse timing data. The implied migra- 
tion rate f into a cylinder with fiducial radius r = 10 10 no cm 
from a disk of outer radius r-i = 10~ 2 r2,2 AU must satisfy 
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Influx rates much larger than these occur in the solar system 
due to a variety of effects that we consider below. 

4. CIRCUMPULSAR DISKS 

A dust disk around the anomalous X-ray pulsar 4U 
0142+61, a magnetar, has recently been detected with Spitzer 
Space Telescope and is estima ted to have mass ~ 10M m and 
a life time of at least 1 Myr iWang. Ch akrabartv & Kaplan 
2006). While other direct evidence for circumpulsar disks 
is lacking, there are several lines of argument that support 
their existence - perhaps transitory — around at least some 
NS. The strongest is the presence of planets around at least 
two NS, one an MSP in the Galactic disk (B 1257+12) and th e 
other in a globular cluster (B1620-26, Sigurdsson et al. 2003). 
For these cases, disks probably arose from ablative disruption 
of a companion star and are not relevant to the model pro- 
posed he re, which concerns cano nical pulsars that are now 
isolated. Michel & Desslerl dl981l) explored how the electro- 
dynamics of pulsars are altered by fossil disks and could alle- 
viate the need for the then "extreme" magnetic field strengths 
of 10 12 G as well as account for some of the observational 
pr operties of pulsars, includ ing nulling. The disks proposed 
by iMichel & Desslerl ( 119811) are substantially more massive 
(up to 10 _4 M Q w 33Mq) than those considered here. 

Fallback disks are plausible features of supernova explo- 
sions because of the wide range of energy and angular mo- 
mentum carried by material that is initially outward-moving. 
Simulations indicate that the fallback disk around a NS typi- 
cally does not provide the conditions needed for planet build- 
ing, primarily because it is born compact and does not spread 
much in radius when it b ecomes neutral and essentially in- 
viscid (Menou et al. 2001). Rocks with large tensile strengths 
can grow from the metal-rich gas and survive well inside the 
R oche radius, r tP (see b elow). 

iMenou et all 120011) and iLu & Chen g (2002) consider 
metal-rich, supernova fallback disks of mass 10 -6 - 10 _3 Mq. 
After a phase of relatively strong accretion onto the NS that 
mimics that in binary systems and produces X-ray emission, 
the disk becomes neutral on a time scale ~ 10 3 — 10 4 yr. The 
disks are born very compact (~ 10 8 - 10 9 cm) and spread to 
10 9 - 1 10 cm by the time th ey become neutral ( Men ou et alJ 
120011) . lEks! & A lpar ( 2005) argue that disks are not disrupted 
by pulsar radia tion pressure and m ay extend into the light 
cylinder. iBlackman & Pernal J2004I) discuss modification of 
the spindown law by an accretion disk and fit models to the 
Crab and Vela pulsars' brak ing indices and s pindown ages to 
constrain disk parameters. Jiang & Li (2005) consider how 
accretion from fallback disks with masses of 10 -6 - 1O _2 M 
will significantly spin down NS during their first 10 4 yr, thus 
providing a mechanism for "injection" of pulsars into the 
Galactic populatatio n over a range of peri ods that extends to 
> 0.1 s ("e.g. lVivekanand & Naravanll98lh . 

Growth of planetesimals in the disk depends on whether 
the material is interior or exterior to the tidal disruption ra- 
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dius, which is, for a self-gravitating object with mass density 
p orbiting a 1.4M NS, 



rtg 



3M* 

27T/9 



1/3 



lO'V'^cm. 



(15) 



The tidal radius is well outside the light cylinder radius, 



r LC =c/n=10 91 Pcm 



(16) 



for all known radio pulsars but is comparable to the light 
cylinders of magnetars with P ~ 10 s. The spin rates of solar- 
system asteroids indicate that t hey are mostly rubble p iles 
with negligible tensile strength ( Richardson et al. 2002) but 
asteroids around pulsars may be stronger macroscopic ob- 
jects if they are metal rich and fractionated and thus can re- 
main intact well inside the gravitational tidal disruption ra- 
dius. A meter-sized rock, R a = 10 2 7?2 cm, with tensile strength 
T = 10 8 7fj dyn cm" 2 , is tidally disrupted at 



GM*pR 2 a 
T 



1/3 



10 



74 cm {pR\T% 1 



,V3 



(17) 



which is about 27 times the NS radius, /?* w 10 km, sim- 
ilar to altitudes w here coherent radio emission is produced 
tovks et al]l2004l) . Thus, should any large ~ 10 km-sized 
planetesimals exist in the disk, they will tidally shred near 
the LCs of 1-s pulsars, while meter-sized rocks can penetrate 
deeply into the magnetosphere until they evaporate. 

5. HEATING OF ASTEROIDS 

Heating of asteroids will ultimately evaporate grains and 
allow them to be ion ized efficiently by radiation from the NS 
and magnetosphere. Cheng ( 1985) concluded that ^0.1 pm 
iron and silicate grains reach distances ~ 1 -3 x lO^cm. De- 
pending on their sizes and inward velocities, meter (or larger) 
sized rocks will reach similar radii before they are evaporated. 

The equilibrium temperature of an asteriod determines 
whether or not it can enter the magnetosphere intact and in 
a largely neutral state. The temperature of a low-albedo aster- 
oid in thermal equilibrium with surface radiation from the NS 



1/2 



: 707 KT 5 r 



■1/2 
10 > 



(18) 



where 7 1 * is the effective radiation temperature that ranges 
from ~ 10 5 to a few times 10 6 K. Cooling mechanis ms are 
not well identified and recent work (Slane et al. 2002) yields 
an upper limit of 10 6 K for the young NS 3C58 while some 
older pulsars show higher temperatures. Moreover, NS mass 
determinations now show a broader range, ex tending from 1.3 
Mq to as high as 2.1 M Q (tNice et al.l l2005). which signifies 
that a broad range of surface temperatures should be expected 
for objects of the same age (Yakovlev & Pethicld 120041) be- 
cause cooling is much more efficient for high-mass NS. 

Heating from magnetospheric radiation is mostly 
by beamed X-rays with luminosity Lx ~ lO -3 ^ 
(Bec ker & Trumped 19971) . 



T aj „^6l2K (g b I 45 P- l5 p- 



i/4 (L x /E 
10" 3 



1/4 



-1/2 



(19) 



where the beaming factor gj, « (S/ir9b) l l A > 1 if the beam of 
size 9b (FWHM) illuminates the asteriod and zero if it does 
not. 



An additional source of heating is ohmic dissipation (in- 
duction heating), discussed in the Appendix. The dominant 
effect is from current flow between the asteroid and magne- 
tosphere. The regime is different from the induction heat- 
ing of aste riods in the solar nebula during th e Sun's T-Tauri 
phase (e.g. Sonettet al. 1970; Herbert 1989), where the am- 
bient plasma density in a much stronger (than present) so- 
lar wind provides a current circuit through the asteroid that 
causes ohmic heating. In the pulsar case, the asteriod is in a 
near vacuum, at least near the light cylinder, so charge carri- 
ers are available primarily from the asteriod itself. In mov- 
ing closer to the NS, evaporation from radiative heating be- 
comes significant (see below), providing a large number of 
charges. During this phase, runaway ohmic heating will occur 
(R. V. E. Lovelace, private communication) even if, as in the 
solar nebula case, ohmic heating is limit ed by diagmagne tism 
from current-induced magnetic fields (Sonett et ali i 19701) . At 
2000 K, which an asteroid reaches at r ~ 10 9 cm around a 
10 5 K and 10 12 G NS, evaporation begins to provide suffi- 
cient charges for significant ohmic heating, leading to run- 
away heating that causes the asteroid to explode. Lower mag- 
netic fields allow the asteroid to reach slightly smaller radii 
before exploding. Perhaps the only way to get asteroids to 
significantly smaller radii is for them to intrude along the 
spin axis of pulsars with nearly aligned magnetic moments. 
The minimum radius reached is perhaps irrelevant as long as 
charged particles are injected inside the magnetosphere where 
they are electromagnetically captured. 

Finally, tidal heating of an object can occur if its or- 
bit is elliptical (or if it inspirals quickly), as occurs with 
sola r-system objects such as I o, Europa and Enceladus (p. 
171 iMurrrav & Dermotilll999l) . iPeale. Cassen & Reynolds I 
( 1979) give the heating rate as 



dE _ 36ire 2 p 2 ( R a V (GM*) 5 / 2 
~d7~~l9~~Qp \r) r l l 2 : 



(20) 



where /i is the rigidity and Q is the specific dissipation rate 
(loss of energy per orbit). The implied radiative equilibrium 
surface temperature is 



10 -1.85^M 
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(21) 



where p\2 = 10~ l2 p. Tidal heating increases faster with de- 
creasing r than other heating processes and could be quite 
large for kilometer-sized objects that reach 10 8 cm; however, 
tidal disruption occurs before this takes place. Thus, unless 
the rigidity is small, tidal heating is at most secondary to the 
other heating mechanims. 

To summarize, radiative heating of asteroids will dominate 
other processes for very hot NS or those with large magne- 
tospheric luminosities beamed toward the asteroid belt. For 
these objects, entry into the magnetosphere will be prevented 
because material is ionized outside the LC. For other ob- 
jects, some combination of radiation and induction heating 
will evaporate asteroids at radii of order 10 9 cm, well inside 
the LCs of many objects. Thus we expect that objects that can 
accrete significant amounts of asteroidal material into their 
LCs will occupy certain regions of the P-P diagram that cor- 
respond to objects with different ages, LC sizes, and ability to 
heat incoming objects. 

For a fixed value of T a at the LC radius, a scaling law for P 
vs. P can be calculated. Using the spindown age, r, = P/2P, 
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as a proxy for chronological age, we find P cx P 5 for magneto- 
spheric or induction heating and P cx p x+l l 2 i for heating from 
thermal radiation whose temperature scales as t s q . Cooling 
curves show a wide range of q, but for objects ~ 10 5 yr old, 
q ~ 0.14 while at ages longer than about 10 6 yr, surface tem- 
peratures f all off steeply, corre sponding to q ~ 0.4 to 1 (e.g. 
Figure 1 of Yakovlev & Pe thick 2004). For a slow falloff with 
age, q » 0. 14, P cx P 46 . Thus heating by magnetospheric radi- 
ation or by thermal NS radiation (in intermediate age objects) 
will yield lines of constant asteroidal temperature at the light 
cylinder distance from the NS that scale roughly as P cx P 5 , 
which are steeper than the putative inner-gap death line. How- 
ever, older pulsars will fall on a locus with P cx P 15 to P 2 3 . 
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FIG. 4. — Asteroid radius, R a , plotted against distance from the NS, r. 
The top horizontal scale is the period for which the light-cylinder radius 
^lc = cP/2-k matches the radius on the bottom horizontal scale and the hor- 
izontal bar designates the range of LC radii for known pulsars. The shaded 
region at the left-hand edge designates the inside of the neutron star. For nom- 
inal parameters, tidal disruption of gravitationally bound objects occurs at 
rtg ~ 10" cm (vertical shaded region). At larger distances, planets may form 
while material closer to the pulsar will only form rocks held together by their 
tensile strength. The leftward going arrow indicates where induction heat- 
ing starts to be important for pulsars with canonical magnetic field strengths. 
Slanted curves to the left of r tg designate rock sizes that can be held together 
by tensile strengths of 10 5 and 10 8 dyn cirT 2 at densities of 8 and 1 gm cnT 3 , 
respectively. The slanted curves to the right are the sizes of planetesimals 
that undergo runaway growth in an annulus of width twice their Hills radius 
in the circumpulsar disk. Vertical dashed lines designate radii at which rocks 
will reach 2000 K for a NS surface temperature of 10 6 K and 2000 K or 3500 
K at 10 5 K. Upward going arrows designate the radii at which an orbiting 
object co-rotates with the NS for spin periods of 0.1, 1 and 10 s. Parame- 
ters include disk surface densities of £ = 1 to 20 gem " 2 (r/10- 2 AUr 3 /2 for 
internal asteroid densities p = 8 to 1 g cnT 3 . 



6. MIGRATION OF DISK MATERIAL INTO THE 
MAGNETOSPHERE 

Disks can supply inward moving asteroids, as we show 
here, that stay neutral until they enter the LCs of cooler and 
longer-period pulsars at rates that are consistent with those 



for nulls and RRATs. While inward migration from radiative 
effects can occur rather slowly for a given object, once it en- 
ters the magnetosphere, we assume that ionized material from 
the asteroids will be subjected to electric forces that drasti- 
cally alters its trajectory. Alternatively, collisions in the disk, 
which occur at interesting rates even for low-mass disks, pro- 
vide more direct injection trajectories into the magnetosphere. 

Studies of accretion onto compact objects usually consider 
ionized gas, which is subject to Lorentz forces that lead to the 
requirement that material needs to orbit faster than the mag- 
netosphere at the Alfven radius in order for it to be accreted. 
Slowl y orbiting material is flung o utward by the propeller ef- 
fect lllilarionov & Sunvaevl |l975l) . resulting in spindown of 
the star. The situation is different if the material is neutral. 
An uncharged rock (or one with small charge-to-mass ratio) 
will not respond to the magnetic field. However, its surface 
charge and polarizability will lead to an interaction. These 
effects are small for macroscopic (meter-sized) objects. The 
main consideration is the set of forces on the rock as it evap- 
orates and becomes ionized. High-melting temperature ma- 
terial can persist to radii of 10 9 1 cm at 2000 K and 10 8 6 cm 
at 3500 K around a NS with surface temperature of 10 5 K 
(Figure |4j. These melting-point radii are similar to those at 
which the rock will co-rotate with the NS, as indicated for 
periods of 0.1, 1 and 10s in Figure |4] By the time the rock 
is ionized, which is probably outside the corotation radius, 

r co = [GM»P 2 /( 27r ) 2 ] 1/3 ~ 10 8 2 cmP 2 / 3 , for most pulsars, it 
encounters regions in the magnetosphere where En jl and is 
electrically captured. 

6.1. Nominal Disk Parameters 

Consider an asteroid disk of total mass M a with minimum 
and maximum radii, r\,r2 r\, and thickness 2H = e/,^ 
yielding a volume Vd w {■K/T)ef l r\, average surface density 

A/ " (22) 



and an optical depth 



3 E 



(23) 



where R a is a characteristic asteroid size and p is its inter- 
nal mass density. For a distribution of asteroid sizes, R a = 
(Rl)/(R 2 a ) which is, for a power-law asteroid size distribution, 
f(R a ) oc Rg" 1 , R a i < Ra < Rai with an equilibrium distribution 
(m = 7/2) (iDohnanvi 111969ft . R a = y/R a \R a 2- The mean time 
between collisions is 

= — , (24) 

where v rms is the rms random velocity. The characteristic 
spreading time for the dis k due to collisions only is (p. 496 
iMurrrav & Dermoti 19991) 



rfi 



(25) 



For nominal parameter values r\ <C ri, r^= 10 V2.-2 AU, 
R a = 1 ^a,km km, and M a = 10" 4 M a ,^M ffi , we have the follow- 
ing typical values for the asteroid disk: 

E = 8 .5 g cm -2 M a - A rf_ 2 (26) 



r=10- 42 (p« a , km )" M a ^rf rl 



(27) 
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f con = 10 32 Gyrp^r^XU^M 

t c = 10 29 yre h pv^ s 4 ,-2^a,km M a,U 

2v z , rms ,„-4 .,,-1/2 3/2 -1 

£/, = ^ = 10 4 U v z ,rm S M^ / 4 r^ rjij 

^V fl = 10 82 M a ,_ 4 (Am)"'- 



(28) 

(29) 
(30) 

(31) 

(32) 



where v rms and v z ,„„ s are rms velocities in km s -1 , with the z 
direction parallel to the disk's axis. We have evaluated f co n 
using an orbital frequency corresponding to an orbital radius, 
r2/2. In these and other calculations we ignore general rel- 
ativistic corrections to orbital parameters, though apsidal ad- 
vance and Lense-Thirring precession will be significant for 
orbits near or inside the LC. However, these effects will likely 
be subdominant to the much stronger electromagnetic forces 
in the region of strong gravity (and strong electromagnetic 
field) close to the surfaces of the NS . 

6.2. Debris Migration 

Collisions can direct some asteroids or their fragments 
into the magnetosphere. However, additional processes now 
known to be significant in the solar system are likely to also 
be important in the circumpulsar disk. These include the di- 
urnal and seasonal Yarkovsky effects for larger rocks (~ 1 
m) and the Poy nting-Robertson effect for small grains (e.g. 
Burn s et alJl979T) . 

The diurnal Yarkovsky effect results from the net force 
on an illuminated object whose afternoon surface is hotter 
than its nighttime surface. It can lead to an increase or de- 
crease in semi-major axis depending on the lag angle of the 
hottest surface. The seasonal Yarkovsky effect applies to ob- 
jects with spin axes tilted from their orbital plane, causing a 
net force along the spin axis that always decreases the semi- 
major orbital radius. In the solar system the diurnal effect is 
larger than the seasonal effect. Poynting-Robertson accelera- 
tion also decreases the orbital radius and is smaller than the 
Yarkovsky effect by a factor v a /c, all else being equal, but is 
important for small grains that are isothermal for which the 
Yarkovsky effect vanishes. Near-Earth objects can migrate up 
to 10~ 3 AU Myr" 1 at 1 AU from the Sun from the seasonal 
Yarkovsky effect, which we take as a fiducial value. 

The Poynting-Robertson effect drags small particles 
into the Sun on a time scale ?pr = Anpc 2 r 2 R/l>L w 
672yrp/?_4r 2 vU L Q /L. NS thermal radiation yields a time 
scale for grains of size lO" 4 /?^ cm at 10 1( Vio cm (e.g. 
iRvbicki & Lightma nll986l p. 153) 



fpR = 16.8 yr 



pR-4r 2 10 



(33) 



For nominal disk parameters, where disk-wide collisions are 
frequent, a steady supply of grains will be injected into the 
magne t osphe re, similar to the interstellar flow proposed by 
lChengllfT985T) . 

We now focus on the seasonal Yarkovsky effect because it 
provides a viable lower bound on the rate of inflow into the 
LCs of some pulsars. The Yarkovsky effect has been detected 
directly through m easured accelerations of asteroids us- 
ing ra dar techniaues iCheslev et all2003tlNesvornv & Bottkd 
2004). It plays a key role in asteroid migration, especially to 
Earth-crossing orbits, by moving objects in the asteroid belt 



into mean-motion resonances that then cause much faster or- 
bital evolution. The seasonal Yarkovsky effect causes a secu- 
lar change in semi-major axis of the form 



roc — 



d-A)L 
cpR^n orh r 2 



AT 

~f~ 



(34) 



where L is the luminosity, p and A are the mass density and 
albedo of the asteriod, and AT /T is the fractional temperature 
difference caused by misalignment of the spin axis from the 
normal to the orbital (or invariant) plane (Rubincam 1998). 
The above proportionality can be derived by assuming sim- 
plistically that a spherical asteroid has a temperature enhance- 
ment AT for the hemisphere facing the direction of motion. 
At r = 1 AU, r w -10~ 3 AU Myr" 1 for a basaltic asteroid of 
100 m size; values of r for iron asteroids are smaller by about 
a factor of two. The effect asymptotes to zero for smaller as- 
teroids as they trend toward isothermality and is smaller for 
more massive asteroids. 

The penetrati on depth for thermal X-rays lx ~ 0.1 — 1/z.m 
dChensIl 9851 IVoil 1 99H) for 0.1 to 1 keV photons, so macro- 
scopic objects will be heated asymmetrically. The ther- 
mal skin depth for iron with thermal conductivity k = 4x 
10 6 ergs _1 cm" 1 K" 1 , heat capacity C p = 5 x 10 6 ergg _1 K" 1 and 
density p = 8 g cm -3 is 



ptt or bC f 
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(35) 



Basalt and other materials have substantially smaller thermal 
conductivities yielding 3 times smaller skin depths, so con- 
ceivable materials can sustain temperature gradients to dis- 
tances near the LC if they are larger than a meter. 

We suppose that rocks possess spins with random orienta- 
tions owing to collisions or chaotic dynamics, so any heat- 
ing will drive the seasonal Yarkovsky effect. In a pulsar disk 
~ 10~ 3 AU from the NS, the Yarkovsky effect will be larger 
than for the solar system for fixed luminosity, but we must 
take into account both the thermal radiation from the NS and 
non-thermal radiation from the magnetosphere. For a thermal 
luminosity Lk, 10 28,9 ergs" 
will be 



1 R 2 * 6^* 4 5' m e Yarkovsky drift rate 
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(36) 



which is much smaller than the reference solar-system rate we 
have adopted, |r ss ,iAu| ~ 10~ 3 AUMyr _1 , owing to the small 
thermal luminosity and the much larger orbital frequency. 

Now we consider the fraction of the pulsar's spindown en- 
ergy loss rate, E, that is in a form that can heat a planetesi- 
mal. The LC is at the transition from the near zone to wave 
zone for rotational energy losses and most of the spindown 
loss rate, E, is in low-frequency waves that will have little ef- 
fect on a rock's dynamics. Cheng (1985) concluded that rel- 
ativistic particles from the pulsar wind heat grains negligibly 
while X-rays dominate the ionization of grains after they en- 
ter the magnetosphere, a result that also applies to the heating 
of grains and rocks in the current model. Typically the non- 
thermal X-ray luminosity of a pulsar satisfies Lx ~ lO -3 ^, so 
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where we have used the beaming factor gt, defined before, 
which can exceed unity if the beam illuminates the asteroid 
but vanishes if not. Thus the Yarkovsky drift rate can be much 
larger for nonthermal radiation than for thermal radiation. 

6.3. Asteroid flux into the light cylinder 

6.3.1. Yarkovsky Migration 

The flux of asteroids through a cylinder centered on the NS 
with fiducial radius r and surface area A c = l-ne^r is N a = 
n a A c \r\ w INa.rfr^ 2 , or 



^ ^±M^r\r[^3-r\r[ 
wlO^'V 1 x 

1/2 2 T 4 
r l0 K *,6 1 *,5 



Kss.lAU 



8. 5 gem" 2 J V10" 3 AUMyr- 



(38) 



(39) 



A rate N a ~ 10~ 3 - 1 s" 1 is probably sufficient to account for 
nulling, so asteroids of 10 m size around NS slightly hotter 
than 10 5 K or a disk more massive than 10~ 4 M ffi can bring the 
rate into this range. 
The lifetime of the disk to inward migration is 



Na 



Td'' 



N a 2r\r\ 

.IQ^Gyr /IQ^AUMyr-n 

about ten times shorter than the collisional lifetime for nom- 
inal parameters. However, for parameter values that make 
the influx rate N a large, the disk lifetime will be substantially 
smaller, but still larger than the characteristic lifetime of a 
canonical radio pulsar. 

6.3.2. Collisional Migration 

A debris disk surrounding a neutron star has many char- 
acteristics similar to a planetary ring system and therefore 
many of the results of ring dynamics can be applied to this 
discussion. Inelastic collisions between particles determine 
the structure and evolution of a disk. On short time scales 
(~ t" 1 orbits, where r is the optical depth of the disk) these 
collisions will result in the asteroids moving towards circular, 
coplanar orbits. On much longer timescales, these collisions 
widen the disk, with most of the mass transported inwards and 
most of the angular m omentum transferred outwards. 

Stewart et al. ( 1984) calculate the inward flux of material 
m in a viscous disk. This hydrodynamical result is consistent 
with a result from the more general approach in which the as- 
teroid behavior is stastically modelled using the Boltzmann 
equation, presented in the same paper. They assumed that the 
asteroids are identical, indestructable, nonspinning spheres, 
which have a size much smaller than the inter-asteroidal spac- 
ing. Furthermore, they assume that the velocity dispersion of 
the particles has a Gaussian form, consistent with simulations 
of ring dynamics. The total flux N a through radius r is then, 



N a - 



9 i i a 

2~i^ a Q~r~dr 



(r 2 fffii/), 



(41) 



where 51 is the Keplerian angular frequency and £ is the sur- 
face density. The viscosity v was calculated numerically by 



iGoldreich & Tremaind dl978l) . by balancing the energy dissi- 
pated during the inelastic collisions and the energy genertion 
due to the viscous stress. 



v w0.15- 
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(42) 



Assuming that S, the velocity dispersion v rms , and the op- 
tical depth t do not vary greatly through the disk, the influx 
rate at the inner edge of the disk (set to be at the light cylinder) 
can be estimated: 

N a = lO^s" 1 -^ r f_ 2 K^M^Ml\ ^v^pi , (43) 

with P in seconds, p in g cm -3 , and v rms in km s" 1 . Substitut- 
ing in the previous expression for optical depth, and assuming 
it is much less than unity, we obtain 



N a = 10 53 *- 1 r£ 2 R:%M;^MZ A p-^ ms Pi 



2., 2 



(44) 



This migration rate is large and would provide material at a 
rate sufficient to disrupt emission processes, even if the mech- 
anism for injecting material into the magnetosphere is ineffi- 
cient. 

7. INSIDE THE MAGNETOSPHERE 

The standard picture is that coherent radio emission in pul- 
sars requires counterstreaming particle flows that include e ± 
pairs produced in cascades. Sustained cascades require accel- 
eration regions in the magnetosphere that yield particles ener- 
getic enough to radiate gamma-rays that then pair annihilate. 
As pulsars spin down, cascades cease as the electric fields in 
acceleration region(s) lessen. However, dormant regions can 
produce triggered pair avalanches that cause coherent radio 
transients. 

7.1. Asteroid Evaporation 

Asteroids rapidly evaporate when they reach radii below 
10 10 cm for r,=5x KTK and below 1(F cm for J* = 10 5 K. 
Figure[5]shows evaporation times for objects of differ ent size 
and e vaporation fluxes assuming carbon composition ICheng l 
Figure[6]shows evaporation times vs. distance from the 
star, assuming that the asteroid surface is in equilibrium with 
surface radiation from blackbody radiation. These figures 
show that evaporation times are comparable to the orbital pe- 
riod only after asteroids are deep within the magnetospheres 
of long-period pulsars. Also, for larger objects (~ 100 m), 
evaporation rates are comparable to the Goldreich- Julian par- 
ticle rate, Nqi, well before they are completely evaporated. 
These results confirm the basic picture that asteroid evapora- 
tion occurs at locations and with rates suitable for modifica- 
tion of magnetospheric electrodynamics. Moreover, the range 
of asteriod sizes expected for a highly collisional asteriod belt 
can account for the broad range of time scales seen in coher- 
ent radio emission. The curves shown in Fig.[6]do not include 
the effects of ohmic heating (mentioned earlier). Ohmic heat- 
ing from v x B currents may be significant (see Appendix) 
during the final stages of evaporation and cause the evapora- 
tion times to vary more steeply with radius than shown in the 
figure. The relative importance of the various heating effects 
will be strongly object dependent, so we defer a detailed com- 
parison to a future analysis. 
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FIG. 5. — Evaporation of material from hot asteroids assuming an evapora- 
tion pr essure for carbon using evaporation pressure data quoted from lChen3 
1 1985). Top: Evaporation times for different sizes. Bottom: Evaporation flux. 
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FIG . 6 . — Evaporation times for asteroids of different sizes for two different 
values of NS surface temperature, as labelled. Also shown (dashed lines) are 
the Keplerian orbital period and the spin period, the latter assuming that r 
corresponds to the light cylinder radius. 



7.2. Electromagnetic Capture 

Ionized evaporant from an inspiraling or infalling asteroid 
is rapidly spun up to the NS's spin rate, f2* owing to the typ- 



ically small electron gyro-radius, r g w 10 2 3 cm (r/R^c) P 3, 
(without relativistic corrections and ignoring asteroidal mo- 
tion). The detailed trajectory of the material depends on its 
field-parallel velocity prior to ionization (e.g. collisional mo- 
tion) and on any parallel electric fields. Because the lifetime 
for many asteroids exceeds a spin period, an asteriod and its 
evaporant will encounter gaps. Ionized material is acceler- 
ated inside such gaps, contributing to the current flow in the 
magnetosphere, some of which is subsequently lost to flows 
outside the LC. The fate of material captured in the closed 
magnetosphere is not clear. It can remain tied to magnetic 
field lines but some may diffuse to larger radii or the open 
field-line region from which it is accelerated to the NS or lost 
by the magnetosphere. 

7.3. Currents and Gaps 

Pulsar magnetospheres will adjust the charge density to «gj 
in order to become force free (i.e. vanishing electric field 
parallel to B, E» = 0). But boundary conditions (or, as for- 
merly assumed, large work functions at the NS surface) pre- 
vent this in "gaps" where En ^ 0. Gap sites include the 
polar cap above the NS surface i n the open field line re- 
gion iRuderman & S utherlandl lT975l) . the slot gap at the in- 
terface between the open and closed field-line regions (e.g. 
I Scharlemann et alii 19781 lAronsl[l983t iMuslimov & Harding I 
2003) and t he outer gap that bounds the null surface fl ■ 
B(r) = (e.g [Cheng et alJl986tlRomani & Yadigaroglu|1995t 
i Cheng et all boOOt iHirotani et al. 1 120031 iZhang et alJ bofrtt 
Takata et alJ2006h . 

Particle acceleration in gaps can yield copious gamma-ray 
emission that drives electron-positron (e ± ) pair cascades if 
there is sufficient residual field En ^ 0, which will depend on 
other sources of charged particles, such as thermionic emis- 
sion from the NS surface or injection of charges from external 
sources. 

We assume that the standard picture is at least approxi- 
mately correct because (a) acceleration of particles to high 
energies in the open field-line region is clearly needed to pro- 
duce beamed radio emission that is consistent with the ob- 
served widths and polarization of pulses; (b) subpulse drift be- 
havior is qualitatively consistent with ExB drifts around the 
magnetic axis, which rely on the existence a gap field; and (c) 
if surface magnetic field estimates are correct, the co nditions 
exist for pair production to occur. Uessner et al.l J200 1 ) discuss 
an alternative to this view. Whether pair production is in the 
form of sustained cascades or in episodic avalanches is open 
to question for particular objects, with the answer dependent 
on NS properties and on external triggers. Since the acceler- 
ating field Em oc p e (r) - /? e ,Gj( r ) m the corotating frame, any 
excess particles injected from outside the LC that diminish 
the charge density difference (or change its sign) must drasti- 
cally alter the population of particles and the fraction of «gj 
needed to do so can be small. Drifting subpulse patterns are 
sometimes coherent, other times disrupted, suggesting that the 
ExB drift rate is time variable. Our model provides a natural 
mechanism for this variability. 

7.4. Coherent Radio Emission 

In our model, neutral particles enter the magnetosphere as 
rocks or grains, dissociate and ionize from thermal and per- 
haps magnetospheric radiation, and perturb the current flow. 
In a pulsar with an active cascade, the intruding charges can 
reduce E\\ and thus the density of cascade e ± pairs. Presum- 
ably this reduces the altitude of radio emission, which is tied 
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to one of the plasma frequencies defined in the system, or 
completely detunes the coherent amplification by altering the 
growth rate of the plasma instability that provides the coher- 
ence. As NS spin down, cascades cannot be sustained and gap 
accelerators become quiescent. Charges from accreted rocks 
may temporarily re-activate the <? ± cascade in one or more of 
the gaps, leading to a burst of coherent radio emission. 

Enough uncertainty exists about the structure of these gaps 
that we sidestep any detailed modeling of the process by 
which deconstructed asteroids disrupt the electrodynamics. 
Disruption undoubtedly depends on NS properties such as the 
surface temperature and spin rate. It also depends on the tra- 
jectories of incoming asteroids, whether they slowly spiral in 
to regions where Eh ^ 0, or alternatively, are on highly eccen- 
tric orbits arising from collisions or resonances outside the 
magnetosphere. 

There appear to be ample configurations and mechanisms 
by which disruption can occur, including (a) direct injection 
of charges into the polar-cap gap and slot gap; (b) perturba- 
tion of the outer gap by itself; and (c) direct injection into 
the outer gap, which triggers an inward moving cascade that 
perturbs the inner gap region(s). We consider the most likely 
case to involve contamination of the outer gap, primarily be- 
cause it borders the LC and its volum e presents a large tar- 
get for incoming material (e.g. Tak ata et alJl200il2006l) . In 
long-period, cool pulsars, the outer gap is unable to sustain 
pair cascades that involve pair production from gap gamma- 
rays interacting with X-rays from the NS surface. However, 
an external flux of particles can trigger an avalanche that, in 
turn, may perturb the inner-gap region of the magnetosphere. 
Such perturbation could enhance the pair production in some 
cases or it could quench the pair production in others. Conse- 
quently, perturbation of the outer gap may lead to a burst or a 
cessation (i.e. a null) of radio emission. 

Perturbation of a gap region by injected charges may in- 
crease pair production in some objects and decrease it in oth- 
ers, depending on whether the region is "starved" of particles 
and on the potential drop across it. In active radio pulsars 
with a steady e ± cascade sustained by a slight deficit from 
the GJ charge density, injected charges will diminish En, thus 
reducing the energies of accelerated particles and the number 
of radiated gamma rays that can pair produce. The pulsar re- 
mains radio loud if the conditions are sustained for driving a 
plasma instability that is responsible for the coherence and if 
coherence occurs at detectable radio frequencies. A change 
in radio emission altitude is expected and it could be lowered 
if it depends solely on plasma densities or it might increase 
if the instability growth rate becomes longer. In objects that 
are mostly radio-inactive because they cannot sustain pair cas- 
cades, injected charges can provide the seed particles needed 
to initiate a short-lived cascade. 

External particles need not provide the full GJ current in 
order to perturb the radio emission, so the actual mass needed 
to produce observable effect s is probably orders of magni- 
tude smaller than given in 33.21 In fact, if the GJ charge 
density were variable by 100%, large torque changes would 
be produced. Depending on the time scales of such variabil- 
ity, torque fluctuations would be manifested as an observable 
contribution to the ubiquitous timing noise found in pulsars. 

7.5. Emission Altitude and Beaming 

Most pulsar models invoke the two-stream instability as 
the source of plasma bunching that underlies the strong co- 
herence of the radio emission (e.g. iRuderman & Sutherland! 



[T^llj^fT987tlAsseo & Melik idze1U998i) . Models differ 
in the nature of the two streams but assume that the instabil- 
ity's growth time must be smaller than the propagation time 
from the NS surface, from which the plasma originates and 
propagates at c. Th is requirement is highly restrictive and 
leads to the proposal (lUsov ll987HAsseo & Melikidze 119981) 
that variability near the NS surface causes particle clouds to 
overtake those from earlier discharges and leads to suitable 
instability growth times. When particles are injected into 
the magnetosphere from asteriods, they provide an additional 
stream of particles that can enlargen the parameter space for 
producing coherent radio emission. 

In the two-stream instability, the maximum growth rate is 
at a frequency proportional to the local plasma frequency, 

1 /2 

which is a function of radius, ui p oc (n± /7I , and which also 
depends on relevant Lorentz factors for the particle streams. 
The equation for a dipolar field line, r{9) = Rsin 2 9/ sin 2 ^), 
where 9q is the angle of a field line at the NS surface, implies 
that emission from relativistic particles is beamed in direc- 
tions parallel to the field line tangents that make an angle 9 /2 
from the radial vector. Emission is therefore beamed into po- 
lar angles 9 e = 39/2, for small angles near the magnetic axis 
and well inside the LC. If the Lorentz factor remains con- 
stant, a change in pair density induces corresponding changes 
in emission altitudes and angles, 

^± = _3^ = - 6 ^ (45) 
n± r 9 e 

where the minus signs signify that an increase in pair density 
induces a reduction in emission altitude and emission angle. 
Asteroidal injection may enhance or diminish the pair density, 
so both increases and decreases can be expected. A given 
object may show a change in pulse profile (mode change), 
become too weak to detect because its beam narrows and thus 
misses the line of sight, or become detectable as a result of a 
wider beam. 

7.6. Subpulse Drift Rates 

Modulation of the charge density by injected charges alters 
not only Eh but also the ExB drift of structures within a gap 
region, such as those of fered as exp lanation s for t he dri fting 
subpulse phenomenon. Ruderman & Suther land! i 19751) de- 
rive a drift velocity for the polar-cap gap region 

n « ^ (f ) V2 « 10 3 W^/^A.^, (46) 

where A$ = 10 12 A$i2 volts is the typical potential drop 
along the field lines in the gap. Injected charges will alter 
this potential drop, so the drift velocity will be affected pro- 
portionally. 

7.7. Spindown Energy Loss 
lChend (fT985) writes the total spindown luminosity as 

Euan = E riA + —3— ULV (47) 

where the first term is the contribution from low-frequency 
dipole radiation and the second term results from the torque 
exerted by magnetospheric currents; /» is the actual current 
and /qj is the Goldreich- Julian current. Cheng considered in- 
terstellar grains that intrude into the LC and contribute to /*. 
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For material originating from a circumpulsar disk, the effects 
are largely the same, but the delivery of extrinsic material and 
its likely episodic behavior are different. Intrinsic contribu- 
tions to /* include ions and e ± pairs. Without an extrinsic 
component, the current may be steady with shot-pulse modu- 
lat ions on the light-travel time across accel erating gaps (<C P) 
or (Cheng & Ruderman 1980; Jones 1982) it may show cycli- 
cal behavior on thermal time scales for the polar cap, which 
Jones (1982) estimates to be several thousand seconds. By 
introducing an extrinsic contribution that has a wide range of 
time scales, the rich variety of radio intensity variations — 
and their variation from object to object — can be explained. 

7.8. Torque Fluctuations and Timing Noise 

As discussed earlier, TOA variations will result from the 
recoil of the NS in response to an ensemble of orbiting aster- 
iods. Accreting asteroids will induce true spin fluctuations of 
the NS through two mechanisms. First, spinup of a mass m a 
through electrical capture in the magnetosphere implies corre- 
sponding spindown of the NS. Second, an increase in current 
from injected charges will produce a torque pulse that acts on 
the NS. We consider these in turn. 

Suppose that an effective fraction f m of the mass m a of an 
asteroid is ultimately spun up into corotation from an orbit 
that is much slower than the NS spin. Once ionized, ions 
and electrons will be accelerated in opposite directions, so the 
spunup mass will be significant only for cases where ions are 
accelerated inward. When ions are accelerated outward and 
lost from the magnetosphere, a torque will be applied to the 
NS if the net acceleration is non-radial, as is likely. Thus, 
we expect in either case that the effective value of f m may be 
sizable. The change in NS spin rate for an electrical capture 
radius, r c , is 

-(n*m a r2/h)(l-{l K /n*), (48) 



where /» is the NS's moment of inertia and fl k is the Keple- 
rian angular velocity, which we assume is much less than f2* 
because objects evaporate well outside the corotation radius. 
The resulting rando m walk in spin rate w ill have a strength 
parameter (e.g. ICordes & Greensteinl 19811) 



S FN = 7e(AlV27r) 2 = n(f m m a r 2 c /h)\tt*/2n) 2 
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where 1Z is the event rate, assumed equal to the asteroid injec- 
tion rate, and parameter values were chosen to yield strength 
parameters and r ms pulse phases similar to those in ferred 
from pulsars (e.g. jCorcles&Helfai^ 
U985tlArzoumanian et alJU994UHobbs et al.ll2004l) . The rms 
pulse phase variation over a time span T is 

v<f>(T) = 
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which is comparable to measured rms timing noise variations 
for nominal parameter values if the effective fraction f m is not 
small. 

The second mechanism, fluctuations in current, 
will alter the torque on a NS to produce spin noise. 



ICordes & Greensteinl J1981I) considered a superposition 
of torque pulses associated with nulling and modeled the 
resulting timing noise, which also would be a random 
walk in the spin frequency. Episodic accretion events 
that produce torque pulses of duration W by a fractional 
change in current, // = <5/*//», yield steps in spin frequency, 
Af2* = fiCiW. The strength parameter for the random walk, 
Sfn = 7ZW 2 ( fiVt 1 '2ir) 2 , corresponds to an rms spin-phase (for 
W = 10 3 W 3 s) 

a^T) » l0- 2 *cy f,(KW) l/2 wl /2 P- l5 p- 2 T y 3 / 2 , (51) 

also comparable to measured timing noise. Asteroidal intru- 
sion into the LC is more likely in older, long-period pulsars 
with small values of P for which a<f, is small, whereas mea- 
sured timing noise is larger for short-period objects with large 
P. It seems likely that there are several processes that con- 
tribute to timing noise, with extrinsic current sources playing 
a major role only in some objects. The recent discovery of 
50% increases in P in pulsar B 193 1+24 jKramer et aUl200fl) 
with W « 10 d, 1ZW ~ 0.2, and // ~ 1 is consistent with this 
picture for timing noise, albeit for a case where the individual 
torque events are easily seen in the timing data. For pulsars 
with faster rates and shorter event durations, W, the torque 
events may be inferred only statistically rather than through 
individual changes in P. 
The ratio of spin-frequency perturbations for the two cases 

is 

_ WhP _ 10 _ 2 fiW 3 h 5 P- 15 



- , (52) 

f m m a r 2 P f m m a A5r 2 lQ P 

Spindown from mass loading dominates for the rather large 
effective asteriod masses assumed (f m m a = 10 f m m a \sg}, 
and the large capture radius; adjustment of these and other 
parameters can modify the ratio by several orders of magni- 
tude either larger or smaller. 

Asteroid-induced torques may also induce changes in spin- 
axis orientation, which would amplify their effects on the 
measured pulse phase. The change in orientation will depend 
on how the instantaneous current density deviates from the av- 
erage current density, so is hard to predict. If the mean torque 
depends on the angle a between the spin and magnetic axes 
as f a , then the induced change in P will be (ICordesll 19931) 
SP ~ S9P(df a /da). Even very small angular changes can 
produce easily detectable, accrued pulse-phase changes from 
this effect. 

8. APPLICATIONS 

8.1. Nulls, Drifts and Mode Changes 

Charges injected from asteroids may be sufficent in number 
that they reduce or shut off the pair cascade and thus alter 
or extinguish the radio emission. Several cases may apply. 
First, injected charges may directly perturb the flow in the 
gap region that is responsible for producing the radio-emitting 
particle flow; it is commonly assumed that this is the polar- 
cap gap at the NS surface. In this case, we would also expect 
the current to increase during a null, leading to an increase in 
spindown torque, as described above. Alternatively, asteroid 
charges are injected into an outer-gap region, producing an 
avalanche of pair production with one charge sign flowing into 
the inner gap. 

All things considered, we expect nulls to be most prominent 
in cool, long-period pulsars for which a wide range of aster- 
oidal conditions can yield influx into their large LCs. Nulls 
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are not impossible in younger, hotter objects, but the require- 
ments are more selective on asteroid size and on trajectories 
into the LC. From a given orbital radius, the LC is a fraction 
~ (cP/Awr) 2 <C 1 of the sphere, so collision fragments large 
enough to avoid evaporation outside the LC must be on highly 
restricted trajectories to hit the target. 

Pulsars that show nulling often also show drifting subpulses 
with variable drift rates and changes in pulse profile (mode 
changes). In some cases, drift rates jump between two or 
more preferred values. We interpret such objects as ones that 
can sustain current flows and pair-production cascades, on av- 
erage, but that are modulated by accreted asteroids according 
to the processes we have identified. Injected charges will alter 
the local difference between the actual and Goldreich- Julian 
charge and current densities, p- pa, which will change the 
parallel electric field, Eh, Boundary conditions at the inter- 
faces between the gap, the closed-field-line region, and the 
NS surface then imply changes in the transverse electric field 
and, hence, in the ExB drift rate. Injected charges there- 
fore can account for variations in drift rate and it is not sur- 
prising to find objects that show both nulls and variable drift 
rates. Why should drift rates have preferred values? This 
is equivalent to asking, why should there be quantized val- 
ues of the parallel electric field? We offer no definitive an- 
swer, but it seems possible that particle source regions in 
the polar-cap region , which drift and themselv es appear fi- 
nite in number (e.g. Desfroande & Rank inl 19991) . may appear- 
and disappear as discrete entities. Alternatively, if the cur- 
rent and particle flows consist of several particle types, with 
one or more depending on injection from external sources, 
the cu rrent and potential drop may show discrete states. Fi- 
nally (Gogoberidze et al. 2005) propose that subpulse drift is 
associated with drift waves in the open field-line region of the 
magnetosphere that modulate the particle accleration. Aster- 
oidal influx, in turn, would modulate this modulation. 

8.2. A Cue-ball Model for Quasi-periodic Bursts and Nulls 
in Bl 93 1+24 

The long ~ 40 d quasi-periodic bursts from B 193 1+24 with 
about 20% duty cycle represent a time scale that is very 
hard to understand from processes internal to the magne- 
tosphere or the NS. Periods of about this duration are as- 
sociat ed with Tkachenk o oscillations of superfluid vortices 
(Coddingt on et al.ll2003l) . but presumably these are cyclical 
with high oscillation "Q" from which it is not clear how low- 
duty cycle bursts would be produced. Orbital processes seem 
richer in possibilities and thus more promising. 

B1931+24 itself is a typical canonical pulsar (P = 0.81, av- 
erage (P) = 10" 141 , and spindown age t s = P/2{P) = \ .6 Myr) 
yet the effects seen from it (bursts and nulls, a 50% increase 
in P during bursts, and no short nulls during a burst) are un- 
common among known objects. Observational incomplete- 
ness due to the impracticality of continuous sampling of large 
regions of the sky may play a role, but the observed effects 
are probably physically rare. Therefore we seek special cir- 
cumstances to explain the observed phenomena. 

8.2.1. Eccentric Asteroid Model 

In the context of our model, the burst period and duty cy- 
cle are naturally explained with a large asteroid of radius 
Rb moving in an eccentric 40-d orbit with semi-major axis 
a = 0.26AU(f O rb/40d) 2/3 . At this distance, R b can be as 
large as 100 km (c.f. Figure |4}. This "cue-ball" or "mini- 



Nemesis" asteroid 3 induces observable effects by deflecting 
smaller asteroids into the loss-cone for light-cylinder insertion 
in much the same way that solar-passing stars perturb comets 
in the Oort cloud. Injection into the magnetosphere would be 
stochastic, so the observed continuity of the bursts, once they 
commenced, suggests efficient injection of a large number of 
objects. Burst durations also would be stochastic. We con- 
sider implausible a model where the cue-ball asteroid directly 
injects charges into the magnetosphere. To do so requires an 
eccentricity very near unity in order that the orbit's minimum 
radius a(l-e) be comparable to the LC radius. Tidal heating 
(c.f. Eq. I20i is also unlikely to be significant for an asteroid 
with 40 d orbital period with reasonable eccentricity. 

During short duty cycle bursts, the spindown rate increases, 
suggesting that the asteroid induces enhancements to the cur- 
rent flow. We associate such episodes to the periastron of a 
large asteroid because, in an elliptical orbit, an object spends 
less time nearer than farther from the star. For example in 
an orbit with e = 0.61, the cue-ball asteroid spends 20% of 
its time at radii less than twice the minimum NS distance 
of a(l-e). If only a 10% increase in distance is sufficient 
to terminate the bursts, an eccentricity of 0.25 is implied. 
The opposite case, where bursts are associated with apastron 
and a termination of asteroidal influx, would imply a lower- 
eccentricity orbit but does not readily explain the increase in 
P during bursts. 

For a large asteriod in the densest part of the disk of small 
asteroids, the mean time between collisions is (c.f. Ea.l30> 

t c =10^s e £^R. d ,rl_ 2 M-U, (53) 

where Av is the relative velocity. This velocity differ- 
ence is, for small asteroids in circular orbits of size equal 
to the periastron distance of the cue-ball asteriod, Av = 
^(1 + e) 1 / 2 - \)(GM/r) ~ (e/2)v, a potentially large velocity. 
For nominal parameters (Eg . I26l32l > t c ~ 1 s for Rb/R a = 10 3 
(i.e. Rb = 1 km), implying a large number of collisions per 
orbit, which may imply that the cue-ball asteriod's orbit will 
evolve too quickly. Thus, either the periastron of the cue-ball 
asteriod is outside the densest part of the asteroid disk or the 
disk is much less dense than our nominal disk. 

The eccentricity of a cue-ball asteriod suggests that the ob- 
ject is part of a family of large objects that co-exists outside 
the small asteriod belt of objects that are injected into the pul- 
sar's magnetosphere. This follows from the need to account 
for the cue-ball's eccentricity. In the solar system, near-Earth 
objects (NEOs) exist as a result of chaotic dynamics in main- 
belt asteroids induc ed by resonances with Jupiter and pos- 
sibly Mars (Murrrav & Dermott 1999). Thus there may be 
additional large objects orbiting B 193 1+24. Given the large 
change in P in the on state for the pulsar, it is possible that 
orbital recoil effects on the TOAs of the pulsar are masked. 
Thus it is of interest to carefully analyze timing data to search 
for the individual effects of multiple orbiting objects. 

8.3. Orbital and Disk Lifetimes 

Requirements on the orbital collision model are that (a) the 
cue-ball orbit must be stable for a long period of time, %, that 

3 Nemesis is the name of the putative solar companion star that was hy- 
pothesized to be the cause for quasi-periodic, terrestrial extinction events 
with period ~ 36 Myr due to perturbation of the Oort cloud of comets i RauD 
11983) . 
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we take to be % = 10 5 fb.s yr. and (b) the influx into the LC 
of the pulsar must be at a rate that can yield continuous and 
sustained bursting over a ~ 8-day period. The first condition 
places a constraint on the change in orbital momentum per 
orbital period based on collisions occuring at a rate f" 1 for a 
fraction of the orbital period, fbPorb.b- This yields 
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The second condition requires that the collisional influx into 
the LC must be at a high-enough rate that the magnetosphere 
receives a continuous flow of evaporated material during a 
burst. Assuming collision fragments are directed isotropi- 

cally, the fraction that enters the LC is /lc ~ (cP/47rr) , 
yielding a mean time between injections of ? c //lc mat we re- 
quire to be less than the evaporation time of an asteroid, f, 
Taken together and using Av(r) /vb(r) ■ 
the collision time using 



evap- 



e/2, we can bracket 
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(55) 



Evaluating for m a /mt, = (R a /Rb) 3 ~ 10~ 15 (e.g. for a 100 km 
cue-ball and R a = 1 m), assuming that collisions take place at 
periastron, r = a(l —e), and using an evaporation time ? e vap = 
10 4 ?evap.4 s (which applies to meter-size rocks for a variety of 
temperatures), we obtain 
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This demonstrates that the orbital configuation can be main- 
tained for a sufficiently long time that we do not require our- 
selves to be observing at a special time. 

8.4. Other possibilities 

Analogs to other solar-system phenomena may apply, 
though they do not readily yield a quasi-periodicity in the 
magnetospheric response. For example, "horseshoe" orbits 
of small asteroids that enclose the L3-L5 Lagrangian points 
in a NS-large-asteroid system have libration periods that can 
be large multiples of the orbital period. These orbits span a 
range of distances from the NS that is wider at the L4.5 points 
than at the L3 point. For orbits that are closer to the NS, one 
might expect an increase in collisions or evaporation rate that 
could enhance injection into the LC. Only if the orbits are 
populated with a small number of objects would we expect to 
see episodic or quasi-periodic effects, however. 

As another variant, the 40 d quasi-period may be associ- 
ated with the period of interaction between two orbiting ob- 
jects in resonance. Mean-motion resonances, such as those 
that produce gaps within the solar system's asteroid belt and 
in planetary ring systems, induce chaotic motions that could 
drive asteriods into the NS. However, this process is likely to 
produce a more continuous and possibly large collision rate 
and, hence, less-intermittent bursts. Consequently, we favor 
the cue-ball model for producing the quasi-periodic effects in 
B 1930+24. 

8.5. Rotating Radio Transients (RRATs) 
RRAT objects identified to date populate the same part of 
the P — P diagram as objects that show prominent nulling, 
with the possible exception of J 1 8 19 — 1458, which has a 
large spindown-derived magnetic field that is similar to that 



of magnetars. We consider RRATs to be associated with 
NS that are similar or identical to canonical pulsars. They 
are observed to be quiescent most of the time either because 
they are intrinsically so or because their radio beams do not 
typically point toward us. RRAT pulses have duty cycles 
that are significan t ly sma ller than those of canonical pulsars 
McLaughl in et alJ (|2006) with comparable periods, suggest- 
ing that the intermittent emission is from a region of the mag- 
netosphere smaller than typical. 

If RRAT objects are active pulsars not usually beamed to- 
ward us, their intermittency could be explained by the injec- 
tion asteroidal material into the outer magnetosphere. The 
pair production rate would then be enhanced, which would 
increase the overall plasma density and the altitude of emis- 
sion. As discussed earlier, this would increase the probability 
of detection and make some objects intermittently detectable. 
Pulses from such objects would be wider than typical if the 
entire open field line region is affected. Only if a restricted 
region were affected might we expect to see narrow pulses. 
This appears to be the case with canonical pulsars for which 
single pulses are almost always narrower than the average of 
many pulses, signifying that the open field line region is not 
instantaneously filled with coherent emitters. 

A perhaps more reasonable alternative is that asteroidal ma- 
terial triggers an inactive outer gap. It is easier to inject ma- 
terial into an outer gap and such gaps should exist in all pul- 
sars. In long period pulsars, however, outer gaps have po- 
tential drops too small to sustain self-generated pair cascades 
( Zhan g et all2004j) . However, external triggers can drive pair 
avalanches in such gaps. It is noteworthy that for objects 
where outer-gap emission is a plausible candidate source of 
radio emission, the radio pulse components are quite narrow. 
The best cases are the Crab pulsar and the MSP B1937+21, 
both of which show high-energy emission that aligns in pulse 
phase (or nearly so) with the radio components. These two 
objects show giant pulses in these same components. Gi- 
ant pulses from the Crab pulsar are very narrow if one cor- 
rects for pulse broadening from radio-wave scattering, e.g. 
W/P ~ 100 us /33 ms ~ 3% The pulses from B 1937+21 
Uenet. An derson & Prince! (12001) are a wider ~ 2.5%. By trig- 
gering the outer gap, pair production and coherent emission 
may be produced directly and we may see the emission from 
either the ingoing or out going beam. In this regard, our model 
is similar to that of Zhang. Gil & Dvks (2006), who postulate 
intermittency and reversal of beaming direction but did not 
propose a specific mechanism for the intermittency, as we do 
here. Alternatively (or in addition), activation of the outer gap 
may induce particle flows that in turn activate or alter pro- 
cesses lower in the magnetosphere, such as the polar cap. 

Our model implies that a triggered avalance would continue 
for only a finite amount of time. However, it can be longer 
than a NS spin period if accretion is sustained. 

In some of the RRATS we may expect to see polar cap emis- 
sion from the pulsar. Thus, a test of the model would be to 
search for this emission. We would expect to see these com- 
ponents at quite different pulse phases. 

9. TESTS AND FURTHER APPLICATIONS 

We have suggested a number of plausbility arguments in 
favor of our model, including the distribution of nulling and 
RRAT pulsars in the P—P diagram (Figure the tempera- 
ture and spin rate regime in which we expect disk material 
can enter a pulsar's magnetosphere in a largely neutral state 
(Figure^, and the distribution of spin-magnetic-moment an- 
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gles in Figure |3 Here we consider ways in which the model 
can be tested and used to infer the properties of material or- 
biting radio pulsars. 

As possible tests, we suggest that identification of counter- 
examples to the trends seen in the P — P diagram would 
be valuable. For example, identifying nulls in MSPs or in 
young pulsars with high thermal luminosity or high mag- 
netospheric luminosity combined with favorable orientation 
would suggest alternative causes in those objects. In the case 
of B 193 1+24, the pulsar showing quasi-periodic bursts with 
correlated increase in P, we suggest that there may be multi- 
ple orbiting objects that may be detectable as arrival time per- 
turbations that must be disentangled from the changes in P. 
A careful investigation of the "timing n oise" in a large sam - 
ple of pulsars, such as that presented in Hobbs et al. (2004), 
may reveal trends in types of timing noise vs. P and P and 
may corroborate the existence of torque fluctuations in those 
objects that also show nulls and other evidence for accretion. 

9.1. Young Pulsars 

Injection of charged particles through asteroid intrusion is 
much less likely from young pulsars that have small magne- 
tospheres (if they have low spin periods), have hot surface 
temperatures (> 10 6 K), and have high non-thermal luminosi- 
ties. Exceptions may exist if occasional, large asteroids on 
low-angular-momentum trajectories can rapidly get inside the 
LC before they are evaporated. Thus, the division between 
"young and hot" and "old and cool" NS may depend on the 
collision rate and size of the asteroid disk, which will deter- 
mine the rate at which events might be seen from a given 
object. Already, the trend appears to be that objects near 
the death line in the P — P diagram are much more prone to 
showing transient effects that can be attributed to asteroids. 
Nonetheless, we may expect to see occasional departures 
from this trend. A notable case may be RRAT J1819- 1458, 
which has a small spindo wn age (1 17 kyr) but has a long pe- 
riod (4.3 s) ( Mc Laughlin etal.l2006l). Thus its LC is large but 
its surface temperature (Reynolds et al. 2006) is also large, 
7; w 0.12 keV w 10 61 K. Larger asteroids (> 100 m) can 
enter, though not deeply, into the magnetosphere before they 
evaporate, allowing the possiblity that they trigger a dormant 
outer-gap region. 

9.2. Millisecond Pulsars 

Most MSPs have white-dwarf companions in compact or- 
bits, so those objects are unlikely to have additional debris 
given their large ages (0.1 - 10 Gyr). Disks may exist around 
isolated MSPs as condensed evaporant from the compan- 
ion star that once provided accreting material responsible for 
spinup of the MSP. Such disks may disappear on time scales 
much shorter than MSP radio-emitting lifetimes. For cases 
where there are disks, MSPs are unlikely to allow injection of 
asteroids into their magnetospheres because their light cylin- 
ders are very small, tlc = 10 74 cm(P/5ms), and the non- 
thermal radiation flux typically is high. 

9.3. Infrared Disk Detection 

Unlike the recent detection of a fairly massive disk around 
the magnetar 4U 0142+61 JWang. Chakraba rtv & Kaplan I 
2006), disks around canonical radio pulsars will be much 
more difficult to detect. 

Brvdenet. al. (2006) recently observed pulsar planetary 
system B 1257+12 in the mid-infrared looking for asteroidal 



dust. No evidence was found for emission, but the limits do 
not rule out a disk comparable to the Solar System's asteroid 
belt. We expect the planetary system around B 1257+12 is less 
likely to contain a massive asteroid belt, because gravitational 
interactions will remove debris at a much faster rate than for 
an asteroid belt without planets. For the objects we consider, 
the disk is much more compact about the central source, so 
the emission should peak in the near infrared. 

We assume that the disk, with mass M a , is composed of 
spherical asteroids of area r a , density p, and albedo A. The 
disk is heated by the thermal radiation of the neutron star (ra- 
dius surface temperature T*). We also assume the disk is 
orthogonal to our line of sight. The temperature T a of the as- 
teroid at distance r from the NS is given by Eq.^D As the disk 
is very small, an observer will see the disk as a point source. 
To calculate the total flux density, we add the flux directed 
towards us from all the asteroids: 



2 dN 
dr7rr a F ufim (T a ) — , 
dr 



(57) 



where dN /dr = 2'nrn a ei l r2 is the number vs. r and we assume 
uniform number density, n a . Adding up the Planck flux from 
each asteroid, we obtain 



45 <7 SB (1-A)7? R\ M a r- dxx 



8tt 5 ^ 2 



•a'2 



1 



(58) 



where x\2 = ( hv/kT*) y/lri^/R*- Using our standard disk pa- 
rameters (Eq. 126132b . the peak flux occurs at a wavelength 

i j 

A ps 13//mr 5 " 1 ^g(r2-2'"i,io) J , when the peak of the integrand 
(x sa 2.8) occurs at the geometric average of the two endpoints. 
The integral is approximately 3 for the standard parameters 
but could be as high as w 6 for asymptotically large disks, 
and will be much smaller for narrower disks. However, nar- 
rower disks would not satisfy the contraints suggested in 33.21 
We find then that the observed flux density is 



(59) 



Though this fiducial value would require long exposures on 
today's best instruments 4 , the uncertainty in the disk parame- 
ters allows a significant increase in fluence. Furthermore, the 
stong sensitivity to the temperature of the neutron star could 
allow for disks to be detected surrounding younger, hotter 
NSs. 

9.4. Radio Disk Detection 

Under suitable conditions, an asteroid belt may be detected 
through reflections of pulsed radio emissi on from asteroids 
larger than the wavelength (Phillipsl ll993l) . The ratio of re- 
flected to pulsed flux is approximately 



S u _ 3M a A 
S p 4irr 2 p 

(Rl 



/A 



Ra- 



Ri 



(60) 
(61) 



where /„ is a geometrical factor that takes into account beam- 
ing of the pulsed flux and the solid angle into which pulsed 
flux is reflected; (RJ) is the average over a distribution of as- 
teriod sizes for sizes larger than the wavelength, implying that 

4 The Spitzer IRAC point-source sensitivity is 0.6 fiiy (la) for a 100 s 
exposure at 3.6 fim i htt p : / /ssc . spitzer . caltech . edu/irac/t 
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the mass M„ in the belt a lso includes only those objects. In his 
analysis, Phillips ( 1 993 ) assumed a fan beam for the pulsar ra- 
diation with f g = l; for many pulsars, a pencil beam appears 
to be more appropriate, leading to f g < 1 . For nominal pa- 
rameters we have 

^ = 1Q-5.2 fAuM a ,-4^ (62) 
Sp P r io^a,km 

For strong pulsars with S p « 1 Jy, the unpulsed reflected flux 
~ 10/iJy. With present instrumentation this is very difficult 
to measure while with the EVLA and the Square Kilometer 
Array, it will be straightforward to detect or to place highly 
constraining limits. 

9.5. High-energy Studies 

Canonical radio pulsars are easily detected in X-rays 
through their magnetospheric emission and, in some cases, 
their thermal surface radiation. For cases where an acceler- 
ation region is reactivated by an asteroid, as is plausibly the 
case for the RRAT objects, we may expect to also see high- 
energy radiation. 

9.6. Interstellar Comets 

While no interstellar comet (ISC) has yet been identified, 
the Oort cloud's existence imp lies that a large number of ISCs 
should exist (e.g. lStern 119901) . A typical density may be n c = 
10 14 «i4pc~ 3 , so for a characteristic velocity of 30 km s -1 and 
a cross section for an asteroid disk of nr\ = \0~ A irr\_ 2 AU 2 , 
the typical collision time is 

f e =10 4 - 6 yr("i4>l_ 2 v3o)~\ (63) 

which is too long to be of interest for our standard disk. 
However, a dense disk extending to r^ ~ 1 AU and having 
high optical depth would show an event about once per year. 
Comet-disk collisions would be brief but might induce highly 
episodic behavior in some objects. 

9.7. Pulsar Dependent Disk Properties 

The accretion model can be used as the basis for inferring 
the presence and properties of disks from the emission prop- 
erties of various objects. Time-variable effects appear in some 
objects and not others, suggesting that there is a wide range 
of disk masses among pulsars and a wide range in instrinsic 
pulsar properties that regulate the inflow of neutral material. 
A variation in disk properties is not surprising because pul- 
sars comprise a runaway stellar population whose velocities 
originate from momentum kicks impart ed at o r shortly after 
the formation of the NS (e.g. lLai et al.ll2001l) . The amount 
of fallback material that ultimately settles into a disk will de- 
pend on the amplitude of the natal kick. A correlation study 
of NS space velocities and their radio intermittency may pro- 
vide additional confirmation of the model and constraints on 
disk properties. 

10. SUMMARY AND CONCLUSIONS 

We have presented a new model that highlights the effects 
that asteroids in circumpulsar disks will have on electrody- 
namic and radiation processes in the magnetospheres of neu- 
tron stars. Implied phenomena may occur in both active radio 



pulsars and those that are usually quiescent. The key element 
is the injection of asteriods as macroscopic objects into pul- 
sar magnetospheres, whereupon they evaporate and provide 
charges that are further processed. The rich set of conse- 
quences of this picture unifies what is otherwise a dizzying 
array of phenomena — drifting subpulses, pulse-shape mode 
changes, nulling, and bursts. Support for the model includes 
the theoretical expectation that orbital debris will accumulate 
in supernova fallback disks and that large asteroids can sur- 
vive for much longer than the typical radio lifetime of a radio 
pulsar. Collisions and radiation-driven migration — processes 
identified in the solar system — inject asteriods into the light 
cylinder and provide charges at rates that perturb acceleration 
regions and presumably radio-emission. Objects that show 
bursts or termination of emission (nulls) with short rise and 
fall times appear to have spin periods and period derivatives 
that are consistent with objects that are expected to be cool 
and have large magnetospheres, just those for which aster- 
oidal contamination is the most likely. Empirical confirmation 
of the model would involve detection of debris disks around 
pulsars that show intermittency and mass-constraining non- 
detections for those that do not. Unfortunately, disk detection 
is unlikely with existing instrumentation and may require the 
EVLA and the SKA to detect reflected pulsar radiation in the 
radio band or a large infrared instrument for direct detection. 
The quasi-periodic burster, B 193 1+24, may provide the best 
opportunity with existing telescopes because the elliptical or- 
bit in our model suggests the presence of additional objects 
that may produce orbital recoil of the pulsar measureable in 
pulse timing data. 

Further studies are warranted to explore asteroidal induced 
torque variations and timing noise, inference of disk prop- 
erties from radio-emission properties, multiwavelength stud- 
ies of objects that do and do-not show associated radio phe- 
nomena, and theoretical study of fallback disks as a function 
of supernova kick amplitudes. Correlation studies will test 
the model and point to other information-gleaning methods, 
such as exploring the relationship between the spin/magnetic- 
moment angle and the occurrence of nulls and finding excep- 
tional objects, such as young, hot neutron stars that show nulls 
and appear to contradict the model. The spin history may de- 
termine which objects retain asteroid belts, as those NS born 
with high magnetic field strengths and high spin rates will sig- 
nificantly ablate asteroids. However, a dense disk may rapidly 
spin down the NS and allow asteroids to form. The durations 
of nulls and bursts are closely related to the lifetimes of as- 
teroids as they enter a pulsar magnetosphere. The lifetime, 
in turn, depends on radiative and particle conditions within 
the magnetosphere, so observable effects may provide robust 
contraints on NS temperatures, radiation mechanisms, parti- 
cle flows, and NS masses via cooling mechanisms that are 
mass dependent. 
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APPENDIX 
INDUCTION HEATING 

We consider two types of induction heating, one from the temporally changing magnetic field experienced by an asteroid that 
has no net current flow through it, the other from a field-induced current through the asteroid. 

Well outside the corotation radius, r co = [GM»P 2 /(27r) 2 ] I//3 pa 10 8 ' 2 cmP 2 / 3 , the asteroid is subjected to a changing magnetic 
field, dB/dt, that is determined by the NS spin; this tends to zero as the asteroid approaches corotation but most asteroids will 
evaporate before reaching corotation. The generated electric field oc dB/dt in the asteroid's frame is dissipated within a skin 
depth 8 <C R a , yielding an equilibrium temperature (with a the conductivity and <tsb the Stefan-Boltzmann constant) 

T a , a =(^-\ l '\pa)-^B^^lOKP-^(Bn^a)^rt 12 (— Y^, (Al) 

V247R7SB/ V^max/ 

where a/a m:ix is the conductivity relative to a maximum that is taken to be that of pure copper, ctcu ~ 5 x 10 17 s" 1 . Asteroidal 
material will have lower conductivities and thus higher temperatures than implied in the equation. However, the skin depth 
S = (c I '2tt) P I fia (~ 1\[P cm for copper) also increases with decreasing conductivity so that an asteroid's resistance reaches 
a maxim um when S ~ 7? a . This limits the heating and thus the temperature coefficient to only a factor of a few larger than in 
Eq. IA1I Ohmic heating increases faster with decreasing radius than radiative heating and would be competitive with radiative 
heating at radii ~ 10 8 cm if asteroids could reach that radius (and ignoring corotation). 

The electric field in the frame of an asteroid moving orbitally or ballistically near the light cylinder of a pulsar with v <C c, 
is E = (SI x r) x B/c, or approximately 10 5 4 volts m~ l Bi2(rLc/r) 2 P~ 3 Rl 6- ^ a circuit is established through the magnetosphere 
and asteroid, the current density J = cE yields ohmic heating power P « R^JE as long as the induced azimuthal magnetic field 
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at the asteroid surface satisfies 5, « 2irR a J /c < B. The implied maximum current density, 7 max ~ cB/2ttR„, corresponds to an 
evaporation flux, 

= *EEL « lO^^s^cm-^nP-^^^c/r) 3 . (A2) 
and a circuit resistance ~ 4r/crLc (cgs) or about 120(r/rLc) Ohms. The maximum heating power is 

^max » c(/? fl B) 2 (r/r LC ) » 10 16 ' 4 erg s" 1 [^ZW^.e] VcA) 5 - (A3) 
When the current density is less than /max, the heating rate is P = (J /J mm )Pmm- The corresponding equilibrium temperature is 

/ P \ 1/4 / P \ 1/4 / P \ 1/4 

-°" K ^(w) (£) ■ 

As an asteroid first enters the magnetosphere at r ~ tlc of a long-period pulsar, it is too cool to provide the necessary particle 
flux (c.f. Figure |5j, so heating is dominated by radiative effects (see text) until evaporation becomes large enough to supply 
the current. For NS with 10 5 K surface temperatures and low magnetic fields, injected asteroids reach ~ 10 9 cm before they 
have 2000 K temperatures (Eq.^]and Figure^, at which their evaporation rates are too low to provide the full current density. 
However, with a canonical field strength of 10 12 G, the asteroid temperature and evaporation rise very rapidly with decreasing 
radius so that runaway heating drives an evaporation avalanche, exploding the asteroid. This occurs at about 10 9 3 cm for a meter- 
sized asteroid. For pulsars with smaller field strengths, this radius will decrease. These numbers are altered if there is significant 
particle flux from elsewhere in the magnetosphere that provides a significant amount of current through the asteroid. 
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Parameters for Nulling and Nulling-Limited Pulsars 
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32 


B92, R93 


J1819+1305 




1.0604 


-15.4 


7.67 


50 




LCX02, ■•• 


J 1946+ 1805 


B 1944+17 


0.4406 


-16.6 


8.46 


50 ±7 




V95, R93 


JUUJ4— U7Z1 


T2fW21 no 
oUUJl— U/ 




1 K A 
— 1 J. 4 


1 .DO 


AA £ J_ 1 1 
44.0 ± l.J 


e: 



vyj, Kyj 


J0754+3231 


B0751+32 


1.4423 


-15.0 


7.33 


34 ± 0.5 


26 


B92, R93 


J1649+2533 




1.0153 


-15.3 


7.46 


30 




L04, ••• 


J0528+2200 


B0525+21 


3.7455 


-13.4 


6.17 


25 ±5 


21 


B92, MR02 


J2321+6024 


B23 19+60 


2.2565 


-14.2 


6.71 


25 ±5 


18 


B92, R93 


J ly+J UU+U 


R 1 QAl nn 


1 A/4 JO 


1 s ^ 


7 AQ 

i .4y 


91 -I- 1 
Zl It 1 


Of, 

zo 


RQO RQ1 

£>yz, Kyj 


J1136+1551 


B1133+16 


1.1880 


-14.4 


6.70 


15.0 ±2.5 


46 


B92, MR02 


J2 113+4644 


B2111+46 


1.0147 


-15.1 


7.35 


12.5 ±2.5 


9 


B92, R93 


J0942-5552 


B0940-55 


0.6644 


-13.6 


5.66 


< 12.5 


25 


B92, R93 


J2330-2005 


B2327-20 


1.6436 


-14.3 


6.75 


12 ± 1 


60 


B92, R90 


juojy+1414 


RnASA_i_1 A 
DU0J0+14 


n iq/iq 


111 
— 1 J. J 


s ns 

J.UJ 


1 -U zl 
1Z it 4 


JU.U 


RQO DQ1 

rsyz, Kyj 


J 1057-5226 


B 1055-52 


0.1971 


-14.2 


5.73 


< n 


90.0 


B92, B93 


J0304+1932 


B0301+19 


1.3876 


-14.9 


7.23 


10 ±5 


38 


B92, MR02 


J2305+3100 


B2303+30 


1.5759 


-14.5 


6.93 


10 ±3 


20.5 


RWR05, R93; 


J 1900-2600 


B 1857-26 


0.6122 


-15.7 


7.68 


10 ± 2.5 


25 


B92, R93 


TOfMR lAIA 
JZU4o~ lOlO 


DZU4J lO 


i .yoio 


1 a n 


0.4J 


1 n n -i- o s 

1U.U It Z.J 


IzL 
J4 


VQS A/IT? no 
v y j, iviKuz 


J0943+1631 


B0940+16 


1.0874 


-16.0 


8.28 


8±3 


53 


B92, R93 


J2157+4017 


B2 154+40 


1.5253 


-14.5 


6.85 


7.5 ± 2.5 


20 


B92, R93 


J0837+0610 


B0834+06 


1.2738 


-14.2 


6.47 


7.1 ±0.1 


30 


B92, MR02 


J0944-1354 


B0942-13 


0.5702 


-16.3 


8.30 


<7 


45 


V95, R93 


T1 ^T)_lT7/1 < 
J 1 JJZ+Z /4j 


D I J JU+Z / 


1 1 O/l Q 

1. lZ4o 


— i j. i 


n 1A 
/.JO 


A J_ o 
it Z 


ZO 


TJOT DOI 

rsyz, Kyj 


J1239+2453 


B 1237+25 


1.3824 


-15.0 


7.36 


6.0 ± 2.5 


53 


B92, MR02 


J0151-0635 


BO 148-06 


1.4647 


-15.4 


7.72 


< 5 


14.5 


B92, R93 


J2022+5154 


B2021+51 


0.5292 


-14.5 


6.44 


<5 


23 


B92, R93 


J 1926+0431 


B 1923+04 


1.0741 


-14.6 


6.84 


< 5 


34.5 


B92, R93 


T1 f- \ AjJYTXI 
J 1014+U / j I 


D lOl z+u / 


1 OflAQ 
1 .ZUOo 


1 /I 
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o.y i 
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RQO DQ1 

oyz, Kyj 


J0953+0755 


B0950+08 


0.2531 


-15.6 


7.24 


< 5 


74.6 


V95,EW01 


J1841+0912 


B 1839+09 


0.3813 


-15.0 


6.74 


< 5 


86.1 


B92, EW01 


J0826+2637 


B0823+26 


0.5307 


-14.8 


6.69 


< 5 


81.1 


B92, EW01 


J1910+0358 


B 1907+03 


2.3303 


-14.3 


6.92 


4.0 ± 0.2 
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B92, R90 
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D 1 Z+u — 0+ 


U.JOOJ 
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— 14. J 


f\ 1 zl 
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J J 


RQO RQ1 

£>yz, Kyj 


J 1456-6843 


B1451-68 


0.2634 


-16.0 


7.63 


< 3.3 


37.0 


B92, R93 


J2317+2149 


B2315+21 


1.4447 


-15.0 


7.34 


3.0 ±0.5 


88 


B92, R93; 


J2022+2854 


B2020+28 


0.3434 


-14.7 


6.46 


< 3 


56 


B92, MR02 


J0152-1637 


B0149-16 


0.8327 


-14.9 


7.01 


<2.5 


84 


V95, R93 


TOO 1 0_i_/IO^/l 

JZZiy+4 /j4 


T300 1 0_i_/IO 
DZ.ll /+4 / 


U. jjoj 


— 14.0 


A /lO 

o.4y 
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A O 

4Z 


T3QO DQ1 

isyz, Kyj 


J08 14+7429 


B0809+74 


1.2922 


-15.8 


8.09 


1.42 ±0.02 


9 


B92, R93 


J0837-4135 


B0835-41 


0.7516 


-14.5 


6.53 


< 1.2 


50 


B92, R93 


J0820-1350 


B0818-13 


1.2381 


-14.7 


6.97 


1.01 ±0.01 


15.5 


B92, R93 


J1932+1059 


B 1929+10 


0.2265 


-14.9 


6.49 


< 1 


35.97 


B92, EW01 


TO 1 1 A 1 /I 
JZ1 10+1414 


RO 1 1 1 /I 
JjZ 1 1 J+14 


n /i /i no 

U.44UZ 


1 S S 

— 1 J. J 


7 is 

/.Jo 


<^ 1 
\ 1 


in n 

JU.U 


RQO DQ1 

rsyz, Kyj 


J 1820-0427 


B 181 8-04 


0.5981 


-14.2 


6.18 


< 0.75 


65 


B92, R93 


J 1948+3540 


B 1946+35 


0.7173 


-14.2 


6.21 


< 0.75 


32.0 


B92, R90 


J1752-2806 


B 1749-28 


0.5626 


-14.1 


6.04 


< 0.75 


42 


B92, R93 


J2055+3630 


B2053+36 


0.2215 


-15.4 


6.98 


<0.7 


34 


B92, R93 


TfMCO 17^Q 

ju4jz— 1 /jy 


oU4JU— 10 


O SzlRQ 


1 zl 
— 14. Z 


0. lo 


\ U.J 


Ozi n 

Z4.U 


RQO DQ1 

oyz, Kyj 


J19 13-0440 


B1911-04 


0.8259 


-14.4 


6.51 


<0.5 


64.0 


B92, R93 


J1823+0550 


B 1821+05 


0.7529 


-15.6 


7.72 


< 0.4 


32 


B92, R93 


J0738-4042 


B0736-40 


0.3749 


-14.8 


6.57 


< 0.4 


17.0 


B92, R93 


J1644-4559 


B1641-45 


0.4551 


-13.7 


5.56 


< 0.4 


33 


B92, R93 


J0630-2834 


B0628-28 


1.2444 


-14.1 


6.44 


<0.3 


13.5 


B92, R93 


J201 8+2839 


B2016+28 


0.5580 


-15.8 


7.78 


< 0.25 


39 


B92, R93 


J1921+2153 


B1919+21 


1.3373 


-14.9 


7.20 


< 0.25 


34 


B92, MR02 


J1534-5334 


B1530-53 


1.3689 


-14.9 


7.18 


< 0.25 




B92, PCE 


J0332+5434 


B0329+54 


0.7145 


-14.7 


6.74 


< 0.25 


32 


B92, MR02 
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Cordes & Shannon 

TABLE A 1 — Continued 



Name (J2000) 


Name (B1950) 


P 


log/> 


logTj 


Nulling Fraction 


OL 


References 






(s) 


(s s _1 1 


(vrl 


<%) 






J 1645-03 17 


B 1642-03 


0.3877 


-14.7 


6.54 


< 0.25 


70 


B92, R93 


TfY7/n 9 SOT 
JU I^L— LoLL 




n 1 a^s 

U. lOOo 


1 1 B 
— LD.O 




\ \}.L 


11 

J 1 


T3QO ~DQ1 


J 1844+ 1454 


B 1842+14 


0.3755 


-14.7 


6.50 


< 0.15 


29 


B92, R90 


J1919+0021 


B1917+00 


1.2723 


-14.1 


6.42 


<0.1 


81 


B92, R93 


J1745-3040 


B 1742-30 


0.3674 


-14.0 


5.74 


<0.1 


51 


B92, R93 


J1731-4744 


B 1727-47 


0.8298 


-12.8 


4.91 


<0.1 


64 


B92, R93 


J lou/— UUJZ 


£5 1 0U*4— uu 


n ao 1 s 


— 1 J. J 


7 1/1 


\ U. 1 


C/ \ 
JU 


T3Q7 A/TP HO 
D7Z, 1V1KUZ 


J0823+0159 


B0820+02 


0.8649 


-16.0 


8.12 


< 0.06 


46 


B92, R93 


J0525+1115 


B0523+11 


0.3544 


-16.1 


7.88 


< 0.06 


78 


B92, R93 


J1935+1616 


B1933+16 


0.3587 


-14.2 


5.98 


< 0.06 


72.0 


B92, R93 


J 1430-6623 


B 1426-66 


0.7854 


-14.6 


6.65 


< 0.05 




B92, PCE 


J0922+0638 


B09 19+06 


0.4306 


-13.9 


5.70 


< 0.05 


48.0 


B92, R93 


J2046+1540 


B2044+15 


1.1383 


-15.7 


8.00 


< 0.04 


40 


B92, R93 


J1559-4438 


B 1556-44 


0.2571 


-15.0 


6.60 


< 0.04 


32 


B92, R93 


J1740+1311 


B1737+13 


0.8031 


-14.8 


6.94 


< 0.02 


41 


B92, R93 


J0629+2415 


B0626+24 


0.4766 


-14.7 


6.58 


< 0.02 


30 


B92, R93 


J0437-4715 




0.0058 


-19.2 


9.20 


< 0.0016 




B92, ••• 


J0835-4510 


B0833-45 


0.089 


-12.9 


4.05 


< 0.0008 


90 


B92, R93 



NOTE. — A list of pulsars known to null and pulsars for which limits have been placed on nulling fraction. The period P, period derviative P, the spin-down 
age T;, and the angle between the spin and magnet ic axe s a are also given. The first re ference i s for t he nulling fraction, and second is a, the angle between 
the ma gnetic and spin axes. R90 refers to lRaiikirl fTWl). B92 to|B iggs 1 1992), R93 t olRankirl HW^\. V95 tolVivekanan3 H^SI). LCX02 tolLorimer et al. I 
1 2006), MR02 tolMitra & Rankinl |2003) . EW01 to lEverett & Weisbergl |200H> . L04 to lLewandowski et alj 120041) . F04 to lFaulkner eTaf] 120041) . and RWR05 
to Redman e t"al] 120051) 



